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We study Brownian motion of a heavy quark in field theory plasma in the AdS/CFT setup 
and discuss the time scales characterizing the interaction between the Brownian particle and 
plasma constituents. In particular, the mean-free-path time is related to the connected 4-point 
function of the random force felt by the Brownian particle. By holographically computing 
the 4-point function and regularizing the IR divergence appearing in the computation, we 
write down a general formula for the mean-free-path time, and apply it to the STU black hole 
which corresponds to plasma charged under three U{1) i?-charges. The result indicates that 
the Brownian particle collides with many plasma constituents simultaneously. 
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1 Introduction 

Brownian motion P-l3] is a window into the microscopic world of nature. The random motion 
exhibited by a small particle suspended on a fluid tells us that the fluid is not a continuum 
but is actually made of constituents of finite size. A mathematical description of Brownian 
motion is given by the Langevin equation, which phenomenologically describes the force 
acting on the Brownian particle as a sum of dissipative and random forces. Both of these 
forces originate from the incessant collisions with the fluid constituents and we can learn 
about the microscopic interaction between the Brownian particle and the fluid constituents if 
we measure these forces very precisely. Brownian motion is a universal phenomenon in finite 
temperature systems and any particle immersed in a fluid at finite temperature undergoes 
Brownian motion; for example, a heavy quark in the quark-gluon plasma also exhibits such 
motion. 

The last several years have seen a considerable success in the application of the AdS/CFT 
correspondence to the study of strongly coupled systems, in particular the quark-gluon 
plasma. The quark-gluon plasma of QCD is believed to be qualitatively very similar to the 
plasma of A/" = 4 super Yang-Mills theory, which is dual to string theory in an AdS black 
hole spacetime. The analysis of scattering amplitudes in the AdS black hole background led 
to the universal viscosity bound [7], which played an important role in understanding the 
physics of the elliptic flow observed at RHIC. On the other hand, the study of the physics 
of trailing strings in the AdS spacetime explained the dissipative and diffusive physics of a 
quark moving through a field theory plasma, such as the diffusion coefficient and transverse 
momentum broadening P4T5]. The relation between the hydrodynamics of the field theory 
plasma and the bulk black hole dynamics was first revealed in [16] (see also [T7]). 

A quark immersed in a quark-gluon plasma exhibits Brownian motion. Therefore, it is 
a natural next step to study Brownian motion using the AdS/CFT correspondence. An 
external quark immersed in a field theory plasma corresponds to a bulk fundamental string 
stretching between the boundary at infinity and the event horizon of the AdS black hole. In 
the finite temperature black hole background, the string undergoes a random motion because 
of the Hawking radiation of the transverse fluctuation modes [T8H20]. This is the bulk dual 
of Brownian motion, as was clarified in [2T|[22]. By studying the random motion of the 
bulk "Brownian string", Refs. [2T1|22] derived the Langevin equation describing the random 
motion of the external quark in the boundary field theory and determined the parameters 
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appearing in the Langevin equation. Other recent work on Brownian motion in AdS/CFT 
includes [23H25|. 

As mentioned above, by closely examining the random force felt by the Brownian particle, 
we can learn about the interaction between the Brownian particle and plasma constituents. 
The main purpose of the current paper is to use the AdS/CFT dictionary to compute the 
correlation functions of the random force felt by the boundary Brownian particle by studying 
the bulk Brownian string. From the random force correlators, we can read off time scales 
characterizing the interaction between the Brownian particle and plasma constituents, such 
as the mean-free-path time tmfp- The computation of tmfp has already been discussed in [21] 
but there it was partly based on dimensional analysis and the current paper attempts to 
complete the computation. 

More specifically, we will compute the 2- and 4-point functions of the random force from 
the bulk and, based on a simple microscopic model, relate them to the mean-free-path time 
tmfp. More precisely, the time scale tmfp is related to the non-Gaussianity of the random 
force statistics. The computation of the 4-point function can be done using the usual GKPW 
rule [316] and holographic renormalization (as reviewed in e.g. [26]) with the Lorentzian 
AdS/CFT prescription of e.g. [271[28]. In the computation, however, we encounter an IR 
divergence. This is because we are expanding the Nambu-Goto action in the transverse 
fluctuation around a static configuration and the expansion breaks down very near the horizon 
where the local temperature becomes of the string scale. We regularize this IR divergence by 
cutting off the geometry near the horizon at the point where the expansion breaks down. For 
the case of a neutral plasma, the resulting mean-free-path time is 

1 

^ Tl^' ^ " ^' ^^-^^ 

where T is the temperature and / is the AdS radius. Because the time elapsed in a single 
event of collision is tcou ~ 1/^? this implies that the Brownian particle is undergoing ~ log A 
collisions simultaneously. (So, the term mean-free-path time is probably a misnomer; it might 
be more appropriate to call t~fp the collision frequency instead.) We write down a formula 
for tinfp for more general cases with background charges. We apply it to the STU black hole 
which corresponds to a plasma that carries three U{1) i?-charges. This is more relevant to 
the actual quark-gluon plasma produced in RHIC and the LHC. 

The organization of the rest of the paper is as follows. In section [21 we start with a 
brief review of Brownian motion in the AdS/CFT setup, from both the boundary and bulk 
viewpoints, taking neutral AdS black holes as simple examples. Then we will discuss Brownian 
motion in more general cases where the background plasma is charged. In section [3l we 
discuss various time scales that characterize the interaction between the Brownian particle and 
plasma constituents. In particular, we introduce the mean-free-path time tmfp, which is the 
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main objective of the current paper, and relate it to the non-Gaussianity of the random force 
statistics using a simple microscopic model. In section HI we use the AdS/CFT correspondence 
to compute the random force correlators that are necessary to obtain tmfp- We present two 
methods to compute the correlation functions. The first one is to treat the worldsheet theory 
as a usual thermal field theory. The second one is to use the standard GKPW prescription 
and holographic renormalization applied to the Lorentzian black hole backgrounds. The 
expressions for the random correlators turn out to be IR divergent. In section O we discuss 
the physical meaning of this IR divergence and propose a way to regularize it by cutting off 
the black hole geometry near the horizon. In section [6l we write down the formula for tmip 
for general black holes and, as an example, compute t^fp for a 3-charge black hole, the STU 
black hole. Section [7] is devoted to discussions. The appendices contain details of the various 
computations in the main text. 



In this section we will briefly review how Brownian motion is realized in the AdS/CFT 
setup [2T|,[22], mostly following [21]. If we put an external quark in a CFT plasma at finite 
temperature, the quark undergoes Brownian motion as it is kicked around by the constituents 
of the plasma. On the bulk side, this external quark corresponds to a fundamental string 
stretching between the boundary and the horizon. This string exhibits random motion due 
to Hawking radiation of its transverse modes, which is the dual of the boundary Brownian 
motion. 

We will explain the central ideas of Brownian motion in AdS/CFT using the simple 
case where the background plasma is neutral. In explicit computations, we consider the 
AdS3/CFT2 example for which exact results are available. Then we will move on to discuss 
more general cases of charged plasmas. 

2.1 Boundary Brownian motion 

Let us begin our discussion of Brownian motion from the boundary side, where an external 
quark immersed in the CFT plasma undergoes random Brownian motion. A general formula- 
tion of non-relativistic Brownian motion is based on the generalized Langevin equation [2^1150] . 
which takes the following form in one spatial dimension: 



where p = mx is the (non-relativistic) momentum of the Brownian particle at position x, and 
' = d/dt. The first term on the right hand side of (12.1 p represents (delayed) friction, which 
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depends linearly on the past trajectory of the particle via the memory kernel '~f{t). The second 
term corresponds to the random force which we assume to have the following average: 

{R{t)) = 0, {R{t)R{t')) = Kit - 1'), (2.2) 

where K{t) is some function. The random force is assumed to be Gaussian; namely, all higher 
cumulants of R vanish. K{t) is an external force that can be added to the system. The sep- 
aration of the force into frictional and random parts on the right hand side of (12. ip is merely 
a phenomenological simplification; microscopically, the two forces have the same origin (col- 
lision with the fluid constituents). As a result of the two competing forces, the Brownian 
particle exhibits thermal random motion. The two functions 7(t) and n{t) completely char- 
acterize the Langevin equation (12. ip . Actually, 7(t) and n{t) are related to each other by the 
fluctuation-dissipation theorem [31] . 

The time evolution of the displacement squared of a Brownian particle obeying (12. ip has 
the following asymptotic behavior |2]: 

{T 
— {t <^ trciax) '■ ballistic regime 
(2.3) 
2Dt (t ^ ^reiax) ■ diffusive regime 

The crossover time scale treiax between two regimes is given by 

1 

trelax = — , 7o = dt-f{t), (2.4) 

7o Jo 
while the diffusion constant D is given by 

D = ^. (2.5) 

7om 

In the ballistic regime, t <^ trciax? the particle moves inertially [s ~ t) with the velocity 
determined by equipartition, |x| ~ y^T/m, while in the diffusive regime, t 3> treiax, the 
particle undergoes a random walk (s ~ \/t). This is because the Brownian particle must be 
hit by a certain number of fluid particles to lose the memory of its initial velocity. The time 
^reiax bctwecu the two regimes is called the relaxation time which characterizes the time scale 
for the Brownian particle to thermalize. 

By Fourier transforming the Langevin equation (12. ip . we obtain 

p{u) = fi{u)[R{u) + K{u)], fi{uj) = . .V ■ (2.6) 

7[u;J — tuj 

The quantity fi{uj) is called the admittance which describes the response of the Brownian 
particle to perturbations. p{u), R{u), K{co) are Fourier transforms, e.g., 

POO 

p{uj)= / dtp{t)e^^\ (2.7) 
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while 7[a;] is the Fourier-Laplace transform: 

poo 

-i[uj] = / dt-f{t)e''^'. (2.8) 
Jo 

In particular, if there is no external force, K = 0, (12. 6p gives 

p{co) = —imujx{u)) = iJ,{uj)R{uj) (2.9) 

and, with the knowledge of /x, we can determine the correlation functions of the random force 
R from those of p or those of the position x. 

In the above, we discussed the Langevin equation in one spatial dimension, but general- 
ization to n = (i — 2 spatial dimensions is straightforwardly] 

2.2 Bulk Brownian motion 

The AdS/CFT correspondence states that string theory in AdS^ is dual to a CFT in (rf — 1) 
dimensions. In particular, the neutral planar AdS-Schwarzschild black hole with metric 



dsl = [-f{r)dt' + {dXr] + ^.dr\ f{r) = 1 - (^)'"' (2.10) 



is dual to a neutral CFT plasma at a temperature equal to the Hawking temperature of the 
black hole. 

In the above, / is the AdS radius, t G M is time, and = (X^, . . . ,X'^~'^) e M"*"^ are the 
spatial coordinates on the boundary. We will set / = 1 henceforth. 

The external quark in CFT corresponds in the bulk to a fundamental string in the black 
hole geometry (12.101) which is attached to the boundary at r = oo and dips into the black 
hole horizon at r = r^', see Figured] The X^ coordinates of the string at r = oo in the bulk 
define the boundary position of the external quark. As we discussed above, such an external 
particle at finite temperature T undergoes Brownian motion. The bulk dual statement is that 
the black hole environment in the bulk excites the modes on the string and, as the result, 
the endpoint of the string at r = oo exhibits a Brownian motion which can be modeled by a 
Langevin equation. 

The string in the bulk does not just describe an external point-like quark in the CFT 
with its position given by the position of the string endpoint at r = oo. The transverse 
fluctuation modes of the bulk string correspond on the CFT side to the degrees of freedom 
that were induced by the injection of the external quark into the plasma. In other words, the 



^We assume that d > 3 and thus n > 1. 
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Figure 1: The bulk dual of a Brownian particle: a fundamental string attached to the bound- 
ary of the AdS space and dipping into the horizon. Because of the Hawking radiation of the 
transverse fluctuation modes on the string, the string endpoint at infinity moves randomly, 
corresponding to the Brownian motion on the boundary. 



quark immersed in the plasma is dressed with a "cloud" of excitations of the plasma and the 
transverse fluctuation modes on the bulk string correspond to the excitation of this cloudj^ 
In a sense, the quark forms a "bound state" with the background plasma and the excitation 
of the transverse fluctuation modes on the bulk string corresponds to excited bound states. 

We study this motion of a string in the probe approximation where we ignore its backreac- 
tion on the background geometry. We also assume that there is no 5-field in the background. 
In the black hole geometry, the transverse fluctuation modes of the string get excited due 
to Hawking radiation [18j. If the string coupling Qs is small, we can ignore the interaction 
between the transverse modes on the string and the thermal gas of closed strings in the bulk 
of the AdS space. This is because the magnitude of Hawking radiation (for both string trans- 
verse modes and the bulk closed strings) is controlled by oc g"^, and the effect of the 
interaction between the transverse modes on the string and the bulk modes is further down 
by gl 

Let the string be along the r direction and consider small fluctuations of it in the transverse 
directions . The action for the string is simply the Nambu-Goto action in the absence of 
a S-field. In the gauge where the world-sheet coordinates are identified with the spacetime 
coordinates x'^ = t,r, the transverse fluctuations become functions of x^: X^ = X\x). 
By expanding the Nambu-Goto action up to quadratic order in X^ , we obtain 



^NC = |rf^xv/-det7,.-^ jdtdr 



{dtX 



I\2 



r'fidrX 



I\2 



= So, (2.12) 



/ 

where 7^,^ is the induced metric. In the second approximate equality we also dropped the 
constant term that does not depend on X^ . This quadratic approximation is valid as long as 
the scalars X' do not fluctuate too far from their equilibrium value (taken here to be X^ = 0). 
This corresponds to taking a non-relativistic limit for the transverse fluctuations. We will be 



^For a recent discussion on this non-Abelian "dressing", see |32) . 
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concerned with the vahdity of this quadratic approximation later. The equation of motion 
derived from fl2.12p is 



[rW + dr{r'fdr)]X' = 0, (2.13) 

where we set X\r,t) oc e"*"^*. Because with different polarizations / are independent and 
equivalent, we will consider only one of them, say X^, and simply call it X henceforth. 

The quadratic action fl2.12p and the equation of motion fl2.13p derived from it are similar 
to those for a Klein-Gordon scalar. Therefore, the quantization of this theory can be done 
just the same way, by expanding X in a basis of solutions to f l2.13p . Because t is an isometry 
direction of the geometry fl2.10p . we can take the frequency u to label the basis of solutions. 
So, let {Ui^{x)}, a; > be a basis of positive- frequency modes. Then we can expand X as 

X\x) = r ^Ku^x) + aluM*]- (2.14) 

Jo ^TT 

If we normalize Mcj(x) by introducing an appropriate norm (see Appendix |A]) , the operators 
a, satisfy the canonical commutation relation 

[a^,a^'] = [al,al,]=0, [a^, a^>] = 27r5{uj - u'). (2.15) 

To determine the basis {uuj{x)}, we need to impose some boundary condition at r = oo. 
The usual boundary condition in Lorentzian AdS/CFT is to require normalizability of the 
modes at r = oo [33j but, in the present case, that would correspond to an infinitely long 
string extending to r = oo, which would mean that the mass of the external quark is infinite 
and there would be no Brownian motion. So, instead, we introduce a UV cut-off§ near the 
boundary to make the mass very large but finite. Specifically, we implement this by means 
of a Neumann boundary condition 

drX = at r = r^-^rH, (2.16) 

where r = Tc is the cut-off surfacejfl The relation between this UV cut-off r = and the mass 
m of the external particle is easily computed from the tension of the string: 



m 



Before imposing a boundary condition, the wave equation (I2.13P in general has two solu- 
tions, which are related to each other hj u ^ —u. Denote these solutions by g±uj{r). They 



^We use the terms "UV" and "IR" with respect to the boundary energy. In this terminology, in the bulk, 
UV means near the boundary and IR means near the horizon. 

^In the AdS/QCD context, one can think of the cut-off being determined by the location of the flavour 
brane, whose purpose again is to introduce dynamical (finite mass) quarks into the field theory. 
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are related by g^jif)* = Q-uii^)- These solutions are easy to obtain in the near horizon region 
r ^ th, where the wave equation reduces to 

{co^ + dl)X^^O. (2.18) 

Here, is the tortoise coordinate defined by 

dr 

= -YTT-y 2.19 



n ~ \, log (^-^ ] (2.20) 



Near the horizon, we have 



{d - l)rH V 

up to an additive numerical constant. Normally this constant is fixed by setting r,,, = at 
r = oo, but we will later find that some other choice is more convenient. From fl2.18p . we see 
that, in the near horizon region r = r//, we have the following outgoing and ingoing solutions: 

guj{r) ~ e*'^''* : outgoing, g-w{r) ^ e~*'^''* : ingoing. (2.21) 

The boundary condition fl2.16p dictates that we take the linear combination 

Ur) = g^r) + e^'-g^^r), e^'^ = -|^^4^. (2.22) 

We can show that 6^ is real using the fact that g_^ = g^. 

The normalized modes u^jit, r) are essentially given by fuiij); namely, u^iit, r) oc e~^'^^ fi^{r) . 
A short analysis of the norm (see Appendix shows that the correctly normalized mode 
expansion is given by 



X{t,r) = r [/.(r)e-^-a. + /.(r)*e^-at] , (2.23) 

th Jo 2tt V2u; 



where fu]{r) behaves near the horizon as 

fu^ir) e''^'^* + r ^ rn (r, -oo). (2.24) 

If we can find such fuj{r), then a,a^ satisfy the canonically normalized commutation relation 

We identify the position x{t) of the boundary Brownian particle with X(t, r) at the cutoff 



x(t) ^ X(t,r,) = r ^^[/^(rje— + /.(r,)*e-*at,]. (2.25) 



The equation fl2.25p relates the correlation functions of xif) to those of a,(^ . Because the 
quantum field X(t, r) is immersed in a black hole background, its modes Hawking radiate [T8] . 

9 



This can be seen from the fact that, near the horizon, the worldsheet action (12.121) is the same 
as that of a Klein-Gordon field near a two-dimensional black hole. The standard quantization 
of fields in curved spacetime [31] shows that the field gets excited at the Hawking temperature. 
At the semiclassical level, the excitation is purely thermal: 

("i"-) - ^^!^- (2.26) 

Using fl2.25p and f l2.26p . one can compute the correlators of x to show that it undergoes 
Brownian motion [21j, having both the ballistic and diffusive regimes. 

In the AdSs {d = 3) case, we can carry out the above procedure very explicitly. In this 
case, the metric (I2.10p becomes the nonrotating BTZ black hole: 

ds^ = -(r2 - rl) de + + dX\ (2.27) 

r"' — rjj 

For the usual BTZ black hole, X is written as X = where = + 2-7?, but here we are 
taking X G M, corresponding to a "planar" black hole. The Hawking temperature (12.111) is, 
in this case, 

T^yf^. (2.28) 
In terms of the tortoise coordinate r^,, the metric (12.270 becomes 

ds' = {r^-r]j){~de + drD+r''dX'', = In f ^1^^ V (2.29) 

The linearly independent solutions to (I2.13P are given by g±uj{r), where 

1 p + tU f p — l\ ' 1 p + IV 



g^{r) = - =1 -— e^"'*- (2-30) 

Here we introduced 

p = — , z/ = — = — . 2.31 
The linear combination that satisfies the Neumann boundary condition (I2.16P is 

drQuirc) l-iv 1 + ipcV f pc-iy (2-32) 



drg-cuirc) l + iu 1-ipciy \pc+l, 

where pc = Tc/th- This has the correct near-horizon behavior (I2.24p too. 

By analyzing the correlators of x{t) using the bulk Brownian motion, one can determine 
the admittance p{uj) defined in (12. 6 p for the dual boundary Brownian motion [21]. Although 
the result for general frequency u is difficult to obtain analytically for general dimensions d, 
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its low-frequency behavior is relatively easy to find; this was done in [21] and the result for 
AdSd/CFTd_i is 



= ^ ^ + 0{y:). 



57r 



(2.33) 



This agrees with the results obtained by drag force computations p 4T0lfT2] . For later use, let 
us also record the low- frequency behavior of the random force correlator obtained in [21] : 



G(«)(ti,t2) = (r[i?(ti)i?(t2)]), 

where T is the time ordering operator. 



IGtt 



(2.34) 
(2.35) 



2.3 Generalizations 

In the above, we considered the simple case of neutral black holes, corresponding to neutral 
plasmas in field theory. More generally, however, we can consider situations where the field 
theory plasmas carry nonvanishing conserved charges. For example, the quark-gluon plasma 
experimentally produced by heavy ion collision has net baryon number. Field theory plasmas 
charged under such global symmetries correspond on the AdS side to black holes charged 
under UiV) gauge fields. 

On the gravity side of the correspondence, we do not just have AdS^ space but also some 
internal manifold on which higher-dimensional string/M theory has been compactified. f/(l) 
gauge fields in the AdS^ space can be coming from (i) form fields in higher dimensions upon 
compactification on the internal manifold, or (ii) the off-diagonal components of the higher 
dimensional metric with one index along the internal manifold. In the former case (i), a 
charged CFT plasma corresponds to a charged black hole, i.e. a Reissner-Nordstrom black 
hole (or a generalization thereof to form fields) in the full spacetime. In this case, the analysis 
in the previous subsections applies almost unmodified, because a fundamental string is not 
charged under such form fields (except for the -B-field which is assumed to vanish in the present 
paper) and its motion is not affected by the existence of those form fields. Namely, the same 
configuration of a string — stretching straight between the AdS boundary and the horizon and 
trivial in the internal directions — is a solution of the Nambu-Goto action. Therefore, as far as 
the fluctuation in the AdS^ directions is concerned, we can forget about the internal directions 
and the analysis in the previous subsections goes through unaltered, except that the metric 
fl2.10p must be replaced by an appropriate AdS black hole metric deformed by the existence 
of charges. 

The latter case (ii), on the other hand, corresponds to having a rotating black hole (Kerr 
black hole) in the full spacetime. A notable example is the STU black hole which is a 
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non-rotating black hole solution of five- dimensional AdS supergravity charged under three 
different U{1) gauge fields [35]. From the point of view of 10- dimensional Type IIB string 
theory in AdS^ x S^, this black hole is a Kerr black hole with three angular momenta in the 

directions This solution can also be obtained by taking the decoupling limit of the 
spinning D3-brane metric [361I38] . Analyzing the motion of a fundamental string in such a 
background spacetime in general requires a 10-dimensional treatment, because the string gets 
affected by the angular momentum of the black hole in the internal directions [TTlEnillD]. So, 
to study the bulk Brownian motion in such situations, we have to find a background solution 
in the full 10-dimensional spacetime and consider fluctuation around that 10-dimensional 
configuration. The background solution is straight in the AdS part as before but can be 
nontrivial in the internal directions due to the drag by the geometry. 

In either case, to study the transverse fluctuation of the string around a background 
configuration, we do not need the full 10- or 11-dimensional metric. For simplicity, let us 
focus on the transverse fluctuation in one of the AdSf^ directions. Then we only need the 
three-dimensional line element along the directions of the background string configuration and 
the direction of the fluctuation. Let us write the three-dimensional line element in general as 

ds^ = -ht{r)f{r)de + ^^dr^ + G{r)dX\ (2.36) 

X is one of the spatial directions in AdS^, parallel to the boundary. It is assumed that 
X(t, r) = is a solution to the Nambu-Goto action in the full (10- or 11-dimensional) space- 
time, and we are interested in the fluctuations around it ^1 The nontrivial effects in the internal 
directions have been incorporated in this metric fl2.36p . We will see how such a line element 
arises in the explicit example of the STU black hole in section [HI In this subsection, we will 
briefly discuss the random motion of a string in general backgrounds using the metric (12. 36 p . 

In the metric (I2.36p . the horizon is at r = th where th is the largest positive solution to 
/(r) = 0. The functions ht{r) and hj.{r) are assumed to be regular and positive in the range 
Th ^ r < oo. Near the horizon r ^ r^, expand /(r) as 

f{r)^2kH{r-rH), A;^ = ^/'(r^,). (2.37) 
The Hawking temperature of the black hole, T, is given by 



(3 27r y hrirn) 
For the metric to asymptote to AdS near the boundary, we have 



1 _ Ih^im) ^2.38) 



/2 ' y ^ 



htf --j^, -j--^ as r ^ oo, (2.39) 



^Note that, under this assumption in a static spacetime, the three-dimensional hue element can be always 
written in the form of (|2.36p . The {t, r) and {t, X) components should vanish by the assumption that X{t, r) = 
is a solution, and the (t, r) component can be eliminated by a coordinate transformation. 
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where we reinstated the AdS radius /. Also, because the X direction f l2.36p is assumed to be 
one of the spatial directions of the AdS^ directions parallel to the boundary, G{r) must go as 

,2 



r 

P 



as 



r — )■ oo. 



(2.40) 



We demand that G{r) be regular and positive in the region vh ^ r < oo. Note that the 
parametrization of the two metric components gu, grr using three functions ht,hr, f is redun- 
dant and thus has some arbitrariness. 

Consider fluctuation around the background configuration X(t, r) = in the static gauge 
where t,r are the worldsheet coordinates. Just as in fl2.12p . the quadratic action obtained by 
expanding the Nambu-Goto action in X is 

1 



5*0 



da'y/^Gg'"'d.Xd^X, 



(2.41) 



where ^^^i, is the t,r part of the metric fl2.36p {i.e., the induced worldsheet metric for the 
background configuration X{t, r) =0), and g = det g^,^. The equation of motion derived from 
the quadratic action f l2.4ip is 




fGX' 



0, 



(2.42) 



where ' = dt, ' = dr- In terms of the tortoise coordinate r^, defined by 

1 [K , 
dr, 



f\ht 

02.421) becomes a Schro dinger- like wave equation 



(2.43) 



X^{r) = 0, 



(2.44) 



where we set X{t,r) = e ^^^ri{r)X^{r) and the "potential" V{r) is given by 



V{r) 



-T] 



dr d 
dr^ dr 



1 dr dr] 
7]"^ dr^: dr 



V 



G-^!\ 



(2.45) 



The potential V^(r) vanishes at the horizon and will become more and more important as we 
move towards the boundary r — oo where Vir) ~ 2r'^ jV^. 

Just as in the previous subsection, let the two solutions to the wave equation fl2.44p be 
g^^if) and g-u^ij) = gujij)* . Near the horizon where V{r) = 0, the wave equation (12.441) takes 
the same form as fl2.18p and therefore 5'±aj('") can be taken to have the following behavior near 
the horizon 



it jojr* 



as r — )■ rn- 



(2.46) 
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If we introduce a UV cutoff at r = as before, the solution fuj{r) satisfying the Neumann 
boundary condition fl2.16p at r = Tc is a hnear combination of g±u]{r) and can be written as 
(12.221) . Using this fui{r), we can expand the fluctuation field X{t,r) as 



= \iWH)['^^7^ [/.(r)e-*a. + /.(r)*e-*at] , (2.47) 



where a^^^a]^ are canonically normalized to satisfy (I2.15p . As before, the value of X(t,r) at 
the UV cutoff r = Tc is interpreted as the position x{t) of the boundary Brownian motion: 
X(t, Tc) = x{t). By assuming that the modes Hawking radiate thermally as in (12.261) . we can 
determine the parameters of the boundary Brownian motion such as the admittance ^{u). 

In general, solving the wave equation (12.441) and obtaining explicit analytic expressions for 
g±u}i fuj is difficult. However, in the low frequency limit — )• 0, it is possible to determine their 
explicit forms as explained in pTJ or in Appendix [B] and, based on that, one can compute the 
low frequency limit of fi{uj) following the procedure explained in [21]. The result is 

fi{u) = -^ + Oiu). (2.48) 

From this, we can derive the low frequency limit of the random force correlator as follows: 

'Gir-H) 



7ra'f3 



+ Oiiu) 



2.49) 



3 Time scales 

3.1 Physics of time scales 

In Eq. (12. 4p . we introduced the relaxation time treiax which characterizes the thermalization 
time of the Brownian particle. From Brownian motion, we can read off other physical time 
scales characterizing the interaction between the Brownian particle and plasma. 

One such time scale, the microscopic (or collision duration) time tcoih is defined to be the 
width of the random force correlator function K,{t). Specifically, let us define 



i-coll 







If K(t) = fi;(0)e~*/*™", the right hand side of this precisely gives tcou- This tcou characterizes 
the time scale over which the random force is correlated, and thus can be interpreted as the 
time elapsed in a single process of scattering. In usual situations, 

^relax ^ ^coU- (3.2) 

Another natural time scale is the mean-free-path time tmfp given by the typical time 
elapsed between two collisions. In the usual kinetic theory, this mean free path time is 
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typically tcou ^ ^mfp ^ ^reiax! however in the case of present interest, this separation no 
longer holds, as we will see. For a schematic explanation of the timescales tcou and tmfp, see 
Figure [2J 

R(t) 

tcoU 

A A 



^mfp 

Figure 2: A sample of the stochastic variable R{t), which consists of many pulses randomly 
distributed. 




3.2 A simple model 

The collision duration time tcoii can be read off from the random force 2-point function K,{t) = 
{R(t)R{0)). To determine the mean-free-path time t^fp, we need higher point functions and 
some microscopic model which relates those higher point functions with tmfp. Here we propose 
a simple model[^ which relates tmfp with certain 4-point functions of the random force R. 

For simplicity, we first consider the case with one spatial dimension. Consider a stochastic 
quantity R{t) whose functional form consists of many pulses randomly distributed. R{t) 
is assumed to be a classical quantity (c- number). Let the form of a single pulse be P(t). 
Furthermore, assume that the pulses come with random signs. If we have k pulses at t = ti 
(z = 1, 2, . . . , fc), then R{t) is given by 

k 

R{t) = Y,^iPit-ti)^ (3-3) 

1=1 

where = ±1 are random signs. 

Let the distribution of pulses obey the Poisson distribution, which is a physically reason- 
able assumption if R is caused by random collisions. This means that the probability that 
there are k pulses in an interval of length r, say [0, r], is given by 

Pk{r) = e-^^^- (3.4) 

Here, n is the number of pulses per unit time. In other words, is the average distance 
between two pulses. We do not assume that the pulses are well separated; namely, we do 

^This is a generalization of the discussion given in Appendix D.l of [2r. For somewhat similar models 
(binary collision models), see [12] and references therein. 
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not assume A ^ If we identify R{t) with the random force in the Langevin equation, 
The 2-point function for R can be written as 



oo 



{Emit')) = E - ^^)^(^' - ^^o).' (3-5) 

k=l ' i,j=l 

where we assumed t, t' G [0, r] and ( )^ is the statistical average when there are k pulses in the 
interval [0,r]. Because k pulses are randomly and independently distributed in the interval 
[0, r] by assumption, this expectation value is computed as 

k 



1 ^ 

= — dti---dtk 

T .In 



. (3.6) 



J2 p{t - u)p{t' - + E ^^^^i)knt - u)p{t' - 1, 

Here, the second term vanishes because {(-i^-j)^ = for i ^ j. Therefore, one readily computes 
V (e,e,P(t - U)P{t' - t,)), = - r dt,P{t - t,)P{t' - t,) 

k 

^- dtiP{t-ti)P{t' -ti). (3.7) 

Here, in going to the second line, we took r to be much larger than the support of -P(t), which 
is always possible because r is arbitrary. Substituting this back into (13. 5p . we find 

/oo 
dtiP{t-ti)P{t' -ti). (3.8) 
-oo 

In a similar way, one can compute the following 4-point function: 
{R{t)R{t')R{t")R{t"')) 



= E ^"''^^ E ^^^^J^^^nP{t - U)P{t' - tj)Pit" - t^)P{t"' - t„)),. (3.9) 

k=l i,j,m,n=l 

Again, the expectation value {eiejem^n)^ vanishes unless some of i,j,m,n are equal. The 
possibilities are i=j^m = n, i = m^j=n, i = n^j = m, and i = j = m = n. Taking 
into account all these possibilities, in the end we have 

{Rm{t')R{nR{n) = {Rmit')R{nRin),,sc + {RmiminRinionn^ (3.10) 

where 

(i?(t)i?(t')i?(t")^r))disc = {Rmit')){R{nR{n) + {Rmin){Rit')R{n) 

+ {Rit)Rit"')){Rit')Rit")), (3.11) 

/oo 
du Pit - u)P{t' - u)P{t" - u)P{t"' - u). (3.12) 
'OO 
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We can think of (13.111) as the "disconnected part" and (I3.12p as the "connected part", or 
non-Gaussianity of the random force statistics. 

In the Fourier space, the expressions for these correlation functions simphfy: 

{R{uJi)R{i02)) = 27r/i(5(wi + i02)P{uJi)P{uj2), (3.13) 
{R{u}i)R{uj2)R{uj3)Ri(^A))disc = (2vryu)^[5(a;i + ^2)5(^3 + ^4) + 5{ui+ uj3)5{uj2 + ^4) 

+ 6{Ui + 004)6(002 + U3)]PicOi)P{u2)PMPM, 

{R{ui)R{u2)RMRM)^^^^ = 27ifiS{ui + U2 + CO3 + uj4)P{ooi)P{oo2)PMP{uji). (3.14) 
In particular, for small Wj, 

{R{uji)R{uj2)) ^ 27ffi5{ui + uj2)P{oo = of (3.15) 

(/?(u;i)i?(a;2)i?(w3)i?(t^4))eonn ~ 27r/i5(u;i + CO2 + CO3 + Ui)P{u = Q)\ (3.16) 

Therefore, from the small frequency behavior of 2-point function and connected 4-point func- 
tion, we can separately read off the mean-free-path time tmfp ~ I/a* and P{(jO = 0), the impact 
per collision. 

The discussion thus far has been focused on the case with one spatial dimension, but 
generalization to n = d — 2 spatial dimensions is straightforward. In this case, the random 
force becomes an n-dimensional vector R^ (t), I = 1,2, . . . ,n. Generalizing (13. 3p . let us model 
the random force to be given by a sum of pulses: 

k 

R\t) = Y,elP{t-U). (3.17) 

i=l 

Here, for each value of i, ef is a stochastic variable taking random values in the (n — 1)- 
dimensional sphere S*"^"^. We also assume that ef for different values of i are independent of 
each other. Then we can readily compute the following statistical average: 

(ejef = — , (^efefef) = r~^^^ • ^^'^^^ 

n n[n + 2) 

From this, we can derive the following i?-correlators: 

{R'iu,)R'iu2)) = ^6''6iui + uo2)PMP{uo2)i (3.19) 
n 

{R\ujr)R\uj2)R''MR''M) = {R\ur)R'{u2)R''MR''M),,,,^ 

+ {R'MR\u2)R''MR''M)^,„ (3.20) 
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where 



{R\uji)R\uj2)R''MR\uj^)) 



disc 



{R\uj{)R\uj2)){R''MR\uj,)) 

+ {R'MR^M){R'MR''M) 

+ {R\uJi)R\uj^)){R\uj2)R''M), (3.21) 



{R\u{)R'{u2)R''{u:,)R\u^)) 



conn 




These are essentially the same as the n = 1 results fl3.13p . f l3.14p and we can compute the 
mean-free-path time tmfp ~ 1//^ from the small u behavior of 2- and 4-point functions. 

We will use these relations to read off tmfp for the Brownian particle in CFT plasma using 
the bulk Brownian motion. 

3.3 Non-Gaussian random force and Langevin equation 

In the above, we argued that the time scale tmfp that characterizes the statistical properties of 
the random force R is related to the nontrivial part (connected part) of the 4-point function 
of R. Namely, it is related to the non-Gaussianity of the random force. Here, let us briefly 
discuss the relation between non-Gaussianity and the non-linear Langevin equation. 

In subsection 12. H we discussed the linear Langevin equation (12. ip for which the friction 
is proportional to the momentum p. In other words, the friction coefficient 7(t) did not 
contain p. Furthermore, the random force R was assumed to be Gaussian. In many real 
systems, Gaussian statistics for the random force gives a good approximation, and the linear 
Langevin equation provides a useful approach to study the systems. However, this idealized 
physical situation does not describe nature in general. For example, even the simplest case 
of a Brownian particle interacting with the molecules of a solvent is rather thought to obey 
a Poissonian than a Gaussian statistics (just like the simple model discussed in subsection 
13. 2p . It is only in the weak collision limit where energy transfer is relatively small compared 
to the energy of the system that the central limit theorem says that the statistics can be 
approximated as Gaussian |13|H^. Furthermore, due to the non-linear fluctuation-dissipation 
relations |15], the non-Gaussianity of random force and the non-linearity of friction are closely 
related. An extension of the phenomenological Langevin equation that incorporates such non- 
linear and non-Gaussian situations is an issue that has not yet been completely settled (for a 
recent discussion, see [B])- 

However, the relation between time scales tcou, ^mfp and R correlators derived in subsection 
l3.2l does not depend on the existence of such an extension of the Langevin equation. Below, we 
will compute R correlators using the AdS / CFT correspondence and derive expressions for the 
time scale tmfp, but that derivation will not depend on the existence of an extended Langevin 
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equation either |3 It would be interesting to use the concrete AdS/CFT setup for Brownian 
motion to investigate the above issue of a non-hnear non-Gaussian Langevin equation. We 
leave it for future research. 

4 Holographic computation of the /^-correlator 

In the last section, can be read off if we know the low-frequency limit of 

the 2- and 4-point functions of the random force. For the connected 4-point function to 
be nonvanishing, we need more than the quadratic term 5*0 in (12.121) or (12.411) . Such terms 
will arise if we keep higher order terms in the expansion of the Nambu-Goto action. This 
amounts to taking into account the relativistic correction to the motion of the "cloud" around 
the quark mentioned in subsection 12.21 In the case of the neutral black holes discussed in 
subsection 12.21 if we keep up to quartic terms (and drop a constant), the action becomes 



where the quadratic (free) part 5*0 is as given before in (I2.12p . 

There are two ways to compute correlation functions in the presence of the quartic term 
(14. 2p . The first one, which is perhaps more intuitive, is to regard the theory with the action 
So + Sint as a field theory of the worldsheet field X at temperature T and compute the X 
correlators using the standard technique of thermal field theory [IH]. The second one, which 
is perhaps more rigorous but technically more involved, is to use the GKPW prescription [HIE] 
and holographic renormalization [26] to compute the correlator for the force acting on the 
boundary Brownian particle. 

The two approaches give essentially the same result in the end, as they should. In the 
following, we will first describe the first approach and then briefly discuss the the second 
approach, relegating the technical details to Appendix |Dl In this section and the next, for 
the simplicity of presentation, we will focus on the neutral black holes of subsection 12.21 

4.1 Thermal field theory on the worldsheet 

The Brownian string we are considering is immersed in a black hole background which has 
temperature T given by (12.111) . Therefore, we can think of the string described by the action 

^More precisely, the computation in subsection l4.2l is independent of the existence of any Langevin equation, 
because we directly compute the R correlators using the fact that the total force F equals R in the m — >■ oo 
limit. On the other hand, in subsection 14. 1[ we compute the R correlators directly, but use the relation l|4.4l) 
derived from the linear Langevin equation. So, the latter computation is assuming that a Langevin equation 
exists at least to the linear order. 




(4.1) 




2 



(4.2) 
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c, 



-T-ip 



C2 



Figure 3: The contour for computing real-time correlators at finite temperature. 



(14.11) just as a field theory of X(t,r) at temperature T, for which the standard thermal 
perturbation theory (see e.g. [46j) is applicable. 

For the thermal field theory described by (14.11) . let us compute the real-time connected 
4-point function 

Gi:Uh,t2,h,h) = (r[a;(ti)a;(t2)x(t3)a;(t4)])eonn 

= (r[X(ti, r,)X(t2, re)X(t3, re)X(t4, re)]),„„,, (4.3) 

where T is the time ordering operator and x{t) = X{t,rc) is the position of the boundary 
Brownian particle. In the absence of external force, K{u) = 0, (12. 6p relates x and random 
force R by 

Riuj) = 4.4 

Therefore, using the low-frequency expression for fi{uj) given in (I2.33p . we can compute the 
4-point function of R from the one for x in (14. 3p . 

As is standard, we can compute such real-time correlators at finite temperature T by 
analytically continuing the time t to a complex time z and performing path integration on 
the complex z plane along the contour C = Ci + C2 + C3, where Cj are oriented intervals 

Ci = [-L,L], C2 = [L,-L], Cs = [-L,-L-t/3] (4.5) 

as shown in Figure [31 L is a large positive number which is sent to infinity at the end of 
computation. We can parametrize the contour C by a real parameter A which increases along 
C as 

Ci : z = X-L (0 < A < 2L) 

C2: z = 3L-X (2L < A < 4L) (4.6) 

C3 : z = -L + i{AL - A) (4L < A < 4L + /3) 

The field X is defined for all values of A. Another convenient parametrization of C is 

Ci: z = t, {-L <t<L), 

C2: z = t, l-L<t< L), (4.7) 

C3 : z = -L-iT, (0<r</3). 
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We will denote by X[j] (i = 1, 2, 3) the field X on the segment Cj parametrized by t and r in 
(14. 7p . Henceforth, we will use the subscript [i] for a quantity associated with Ci. 

The path integral is now performed over X[2](t), and X[3](r), but in the L — )■ oo 

limit the path integral over X[3] factorizes and can be dropped |12j- Therefore, with the 
parametrization (14. 7p . the path integral becomes 

VX e'^ j I?X[2] e*(^iii-^Pi), (4.8) 

where S'[j], 2 = 1,2 are obtained by replacing X with X[j] in f l4.ip . The negative sign in front 
of S'[2] in fl4.8p is because the direction of the parameter t we took in f l4.7p is opposite to that 

of C2. 

The correlator fl4.3D can be written as 



G^^Atu t2, h, U) = (rc[X[i] (ti, re)X[i] (t2, re)X[i] (t3, r,)X[i](t4, rc)]),„,„, (4.9) 

where Tc is ordering along C (in other words, with respect to the parameter A), and can be 
computed in perturbation theory by treating Sq as the free part and S'int as an interaction. 
In doing that, we have to take into account both the type-1 fields and the type-2 fields 
X\2]- Namely, we have to introduce propagators not just for but also between and 
X[2] as follows 

D[n](t-t',r,r') = {Tc[X^iit,r)X^^0 ,r')]) ^ = {T\X{t,r)X{t' .r')]), = Dpit - t' ,r,r'), 

D[2i] (t - t', r, r') = {Tc[X[2] (t, r)X[i] (f, r')])^ = (X(t, r)X(t', r')]), = Dw{t - t' , r, r'). 

(4.10) 

Here, ( )o is the expectation value for the free theory with action 5*0 at temperature T. We 
see that the propagators D[ii] and D[2i] are equal, respectively, to the usual time-ordered 
(Feynman) propagator Dp and the Wightman propagator D\y of the field X{t,r). We must 
also remember that we have not only interaction vertices that come from S^^^ and involve 
X[i], but also ones that come from S^2] involve X[2]. The second type of vertices come 
with an extra minus sign. 

Using the propagators fl4.10p . the connected 4-point function is evaluated, at leading order 
in perturbation theory, to be 



G'conn('^b ^2, ^3, ^4) = y^2t^5{uji + ■ ■ ■ + J dr 



perm 

(ijki) 



^'^^ {uJi)Din] (wj) + r'^fdrDiu] {uji)drD[n] {^^j. 



f 

f 



A21]) • (4.11: 
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Here, we wrote down the result in the Fourier space and used a shorthand notation D[ii]{ui) = 
D[ii]{uJi,r,rc). The summation is over permutations (ijkl) of (1234). 

We are interested in the low frequency limit of this correlator. In that limit, the propaga- 
tors simplify and can be explicitly written down. In Appendix [CI we study the low-frequency 
propagators, and the resulting expressions are 



(4.12) 



r| cc;(l - e-/3-) ' 

where r* is the tortoise coordinate introduced in fl2.19p . As explained in (IB.2ip . the precise 



low frequency limit we are taking is 

uji — )■ 0, /3, cjjn : fixed. (4-13) 

The reason why we have to keep Uir^, fixed is that, no matter how small Ui is, we can consider a 
region very close to the horizon (r^, = — oo) such that cjjr^, = 0{1). If we insert the expressions 
fl4.12p into (14. lip and keep the leading term in the small Ui expansion in the sense of (I4.13p . 
we obtain 

<^conn(^l> ^2, ^3, W4) ^^^^^ 5{Ui H h ^4) 

X ^ {uji + u,) C dr/je-^"^'^^+''^'>'* + 0{u-^), (4.14) 

l<i<j<4 -^^"^ ^ 

where we ignored numerical factors. Using (14. 4p and (I2.33p . we can finally derive the expres- 
sion for the R correlator: 

X ^ (c^, + c^^.) r rfr,ye-2'("'+">*+0(a;2^. (4.15) 

l<i<j<A ~°° ■' 

Let us look at the IR part of (I4.14p . namely the contribution from the near-horizon region 
(large negative r*). Because / ~ (ci— X'^e'^'^^^^'^"^* near the horizon, the r^, integral in (I4.14p is 

oo J d 1 J_oo 

which diverges because of the contribution from the near horizon region, r* — t- — oo. We will 
discuss the nature of this IR divergence later. 



22 



4.2 Holographic renormalization 

Next, let us discuss another way to compute the correlators of the boundary Brownian motion, 
following the standard GKPW procedure [5l|6]. For this approach, we send the UV cutoff 
— )■ oo and let the string extend all the way to the AdS boundary r = oo. The boundary 
value of X{t, r) is the position of the boundary Brownian particle: x{t) = X{t, r — t- oo). The 
boundary operator dual to the bulk field X{t,r) is F(t), the total force (friction plus random 
force) acting on the boundary Brownian particle. The AdS/CFT dictionary 

\ I CFT 

says that, to compute boundary correlators for F, we should consider bulk configurations for 
which X(t, r) asymptotes to a given function x(t) at r = cxd, evaluate the bulk action, and 
functionally differentiate the result with respect to x(t). Note that, in the limit Tc — )■ oo or 
m — > oo that we take, friction is ignorable as compared to random force i?, and F correlators 
are the same as R correlators |12j . Roughly speaking, because the Brownian particle does not 
move in the m — > oo limit, there will be no friction and thus R = F. 

In the end, the resulting 4-point function (FFFF) is essentially given by the interaction 
term in the action, with the X fields replaced by the boundary-bulk propagators. Namely, 



{r[F{h)F{h)F{t,)F{U)]) ~ / dtdr 



perm 

(ijkl) 



X 



/ 

dtK{tk) dtKit 



+ rV drKiti) drK{tj 



+ r^fdrK{tk)drK{ti 



(4.18) 



/ 

where K{ti) = K{t,r\ti) is the boundary-bulk propagator from the boundary point tj to the 
bulk point {t,r). This is the Witten diagram rule that we naively expect. However, because 
the worldvolume theory of a string is different from, e.g. a Klein-Gordon scalar, a careful 
consideration of holographic renormalization (see |26] for a review) is necessary. Indeed, 
the naive expression is fl4.18p is UV divergent and needs regularization. Furthermore, our 
black hole spacetime is a Lorentzian geometry and we should apply the rules of Lorentzian 
AdS/CFT [STJEHj. As is explained in Appendix [D| after all the dust has been settled, the F 
correlator gives exactly the same IR divergence as the naive computation of the R correlator, 
(14.151) . This implies that this IR divergence we are finding is not an artifact but a real thing 
to be interpreted physicallylfl 

It is worth pointing out that the result (14.181) has a similar structure to the one we saw in 
the toy model (13.121) . with the propagator K{t) roughly corresponding to P{t). It would be 
interesting to find an improved toy model which precisely reproduces the structure (I4.18p . 



^Although the IR parts are the same, the resuh obtained in the previous subsection 14 . 1 1 using the worldsheet 
thermal field theory is not quite the same as the one obtained in this subsection 14.21 using holographic 
renormalization, due to the counter terms added to the latter at the UV cutoff r ~ Vc- 
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In Appendix ID.5t we also computed the retarded 4-point function of random force. The 
expression is free from both IR and UV divergences and the final result is finite. However, 
because we do not know how to relate the retarded 4-point function and tmfp? this cannot be 
used to compute tmfp- It would be interesting to find a microscopic model that directly relates 
retarded correlators and tmfp- 



4.3 General polarizations 

The argument so far has been as if there were only one field X and the associated random 
force R. However, in the general d > 3 case we have n = d — 2 > 1 fields , / = 1, 2, . . . , n. 
Considering all X^ , the bulk action ( 14. 2 p actually becomes 

1 2 



5i 



int 



16na' 



dt dr 



I\2 



f 



Il\2 



(4.19) 



The associated random force has n components too. 

The computation of 4-point functions in this multi-component case can be done completely 
in parallel with the one-component case. Let us define 



G^:}J^''\h,h, h, U) ^ {T[X\U, r,)X'\h, r,)X^(t3, r,)X^(t4, r,)]). 



(4.20) 



This is nonvanishing only if some indices are identical. More precisely, the only nonvanishing 
cases are (i) all indices are identical, I = J = K = L, or (ii) indices are pairwise identical, 
I = J^K = L, I = K^J = L, or I = L^J = K. 

In case (i), the resulting 4-point function is exactly the same as the one-component case 
(14.111) . Consequently, the IR form of the random force correlator Gconn^^^ is the same as the 
one-component case (14.151) . 

In case (ii), on the other hand, the 4-point function becomes 



UJIUJ2 



-2T[5{ijJi 



Ui) I dr 



X 



/ 



IGvra' 

Dfii] {ui)Dui] M + r^fdr-Diu] iuJi)drD[u] (^2) 



X 



/ 



Diii^{uj3)Diu]{uJi) + r^/<9,.L'rii](a;3)(9,.L'rii](a;4) 



- (-D[ii] D121], 



(4.21) 



The IR form of the random force correlator is 



Gioli^'i^U C^3, C^4) ~ ^^(^1 + ■ ■ ■ + W4) 



X ^ {uji + ujj)e 

l<j<2,3<j<4 



^2 

dr^ — 

00 J 



(4.22) 
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Comparing this with the expectation from the field theory side, (13.221) we observe the same 
structure. Namely, the connected 4-point functions are nonvanishing only when the polar- 
ization indices are all or pairwise identical. The precise relative values of the nonvanishing 
4-point functions are model-dependent and not important; in the simple model of subsection 
13.21 it depends on our choice of the expectation values (13.181) . 



5 The IR divergence 

In the last section, we computed the connected 4-point function for the random force R and 
found that the low-frequency expression, 

G'conn(^l> ^2,^3,^4) ~ ^^5(u;i + --- + W4) 



X 

i<j<i<4 

has an IR divergence coming from the integral in the near horizon region. What is the physical 
reason for this divergence? Very near the horizon, the expansion of the Nambu-Goto action 
in the transverse fiuctuation X breaks down because the proper temperature becomes higher 
and higher as one approaches the horizon and, as a result, the string fiuctuation gets wilder 
and wilder. The correct thing to do in principle is to consider the full non-linear Nambu-Goto 
action, but this is technically very difficult. Instead, a physically reasonable estimate of the 
result is the following. Let us introduce an IR cutoff near the horizon at 

rs = rH + e, (5.2) 

where e <^rH- We take this cutoff to be the radius where the expansion of the Nambu-Goto 
action becomes bad. Then, in IR-divergent expressions such as (14.151) . we simply throw away 
the contribution from the region r < r < by taking the integral to be only over r > r^. Of 
course, to obtain a more precise result, we should include the contribution from the region 
th < r < Vg with the higher order terms in the expansion of the Nambu-Goto action taken 
into account. However, we expect that the contribution from this region th < r < will be 
of the same order as the contribution from the region r > and, therefore, we can estimate 
the full result by just keeping the latter contribution. 

With this physical expectation in mind, let us evaluate the mean-free-path time tmtp by 
introducing the IR cutoff (15. 2p . The parameter e appearing in (15. 2 p can be related to the 
proper distance from the horizon, s, as follows: 



(5.3) 



v7 Jrn ^{d-l)rH{r-rH) \{d-l)rH 
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Therefore 



e ~ s'rn, (5.4) 
where we dropped numerical factors. In the tortoise coordinate r^, the cutoff is at 

\. logg', (5.5) 

where we used (12. 20 p . 

The introduction of an IR cutoff of the geometry near the horizon also means that the 
resulting expressions such as (15. ip . with the IR cutoff imposed, is valid only for frequencies 
larger than a certain cutoff frequency Ug. We can relate Ug with the geometric cutoff as 
follows. If we cut off the geometry at r^, = r^, we have to impose some boundary condition 
there (just as for the brick wall model). For example, let us impose a Neumann boundary 
condition. As was shown in (IB.lQp . for very low frequencies, the solutions to the wave equation 



behave as 

/.(r)~e*"'-*+e-*-'-*. (5.6) 

For this to satisfy Neumann boundary condition dr,fuj{r)\r,=ri = 0, we need u = nir/rl where 
n E 1j. Namely, the frequency has been discretized in units of vr/|r^|. Therefore, the smallest 
possible frequency is 

|r|| f3log{l/s) 

If we use fl5.5l) and (15.71) . the correlator fl5.ip becomes 
Gi±{oo,,oo2, cos, io,) ~ -^6{oo^ + ■ ■ ■ + u,) oo^rl [ dne-^'^'>"'- 



7^'^('^i \-UJi)uJsrH 

■ 1 ■ 1 

~ ^^(^1 + ■ ■ ■ + C.4) ~ ^,^5 1og(l/,)^(^^ + ■ ■ ■ + C.4). (5.8) 

On the other hand, from fl2.35p . the 2-point function is 

G^''\uji,uj2) ~ -^6{uJi + uj2) (5.9) 

Comparing above results and the toy model results fl3.15p . f l3.16p . we obtain 

1 a'P 



— 



P(w = 0) , (5.10) 



^i s2 1og(l/s)' l3s^\og{l/s)' 
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Now the question is how to determine the length s. This must be the place where the 
expansion (14.11) of the Nambu-Goto action becomes bad. One can show that this occurs a 
proper length ~ y/a' away from the horizon due to thermal fluctuation (Hawking radiation) 
in the black hole background (for an argument in more general setups see subsection 16. ip . 
This leads us to set 

s-Vcy. (5.11) 

At this point, the local proper temperature becomes of the order of the Hagedorn temperature, 
~ l/y/a'. The above condition must be the same as the condition that the loop correction 
of the worldsheet theory to the 4-point function {F^) becomes of the same order as the tree 
level contribution. 

If we substitute (15.111) into (15.101) . we obtain 

1 TA^/^ 
t^fp ~ — — , P{uj = 0) ~ -y== (5.12) 
I log A V logA 

where, following the convention of the d = 5 (AdSs) case, we defined the '"t Hooft coupling" 
by 

A (5-13) 
a' 

where we restored the AdS radius / which we have been setting to one. 

The result (I5.12p is quite interesting. In [2T], the collision duration time tcou was deter- 
mined to be 

tcoii ~ ^. (5.14) 

Therefore, tmfp given in (I5.12p implies that a plasma particle can be thought of as in interaction 
with roughly log A other particles simultaneously. 

Even if we take into account the fact that has in general more than one component 
(/ = 1,2, ... ,n = d — 2) and use the results such as (I3.22p . (14.221) . we end up the same 
estimate for t^fp as far as its order is concerned. 

6 Generalizations 

In the previous section, we derived using AdS/CFT the expression for the mean- free-path 
time tjnfp for the simple case of neutral plasma. In this section, we sketch how this generalizes 
to the more general metric (I2.36P and present the expression for the mean-free-path time for 
more general systems such as charged plasmas. As an example, we will apply the result to 
the STU black hole. 
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6.1 Mean- free-path time for the general case 

We are interested in computing the mean-free-path time in field theory by analyzing the 
motion of a Brownian string in the metric fl2.36p . For that, as has been explained in section [3] 
for the neutral case, we need to compute the 4-point function of the random force in addition 
to the 2-point function. 

Expanding the Nambu-Goto action in the background metric fl2.36p up to quartic order, 
the action for the string in the tortoise coordinate defined in fl2.43p is given as follows: 



5" = 5*0 + 5'int, (6.1) 

So = -^, [ dtdnGiX^-X'^), (6.2) 
Area J 

1 /■ . . 



Sint = I dt dr, f-(X2 - X'2)2, (6.3) 

IGna' J htf 

where we dropped a constant independent of the field X, and ' = dt, ' = dr^. As we discussed 



in subsection 14.21 for the simple neutral case, we can use S'int as the interaction term and 
apply the usual GKPW rule to compute correlators for the random forcelf F dual to the bulk 
field X. As before, the naive result from the GKPW prescription includes both UV and IR 
divergences. Using holographic renormalization, which is discussed in Appendix [D] for the 
neutral case, we can remove the UV divergence by adding counter terms to the action. The 
IR divergence, on the other hand, signals the breakdown of the quartic approximation f l6.ip . 
We regulate this divergence by introducing an IR cutoff at r^, = near to the horizon, whose 
physical motivation was explained in section [51 

Following the same analysis as in section |5] now with the interaction term (16. 3p . we obtain 
an expression similar to fl4.15p for the connected random force 4-point function. The dominant 
contribution comes from the near-horizon region and is given in frequency space by 

(r[Fl)con„ ~ + ■ ■ ■ + ^4) / dn^ Yl + ^.)e-^^(--^-^)^% (6.4) 

^ •^I'l J ^ l<i<j<A 

where is the aforementioned IR cutoff (in the tortoise coordinate). Let the IR cutoff in the 
r coordinate be at r = r/^ + e = rg. The parameter e is related to the proper distance s from 
the horizon as 

2thr{rH) s'^kn , . 

Using the relation fl2.43p between and r^, we can estimate the cut-off integral (16. 4p as 

{T[F^]) - ^^^^5(a;i + --- + uj,), (6.6) 




^Recall that in this setup the force F is equal to the random force R. 
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where Ug is the smallest frequency for which the expansion (16. ip is valid. Combining this with 
the result (12.491) for the 2-point function, the mean-free-path time is estimated as 

trafp ~ 5 • (6.7) 

Now, let us determine the IR cutoff parameters s (or equivalently e) and Ug appearing in 
(16.71) . As before, we take the IR cutoff to be the location where So and Si^t become of the 
same order. As is clear from (16.21) . (16. 3p . the expansion of the Nambu-Goto action becomes 
bad at the location where 

• ^^-^ - ' 

So, we would like to estimate X, X'. Near the horizon, r ^ r^, we can write the action (16. 2p 
as 



Sor^^l dtdnix' -X"), X^^^^X. (6.9) 

There being no dimensional quantity in the problem other than the temperature T, we must 
have X,X' ~ T, namely \X'\ ~ /G{rH) T. So, the condition (16. 8p determines the 

IR cutoff to be at 

r-rH = e~— — — -. 6.10 

kuhtirH) 

In term of s, the IR cutoff is at the string length: 

s-Vo'. (6.11) 

It is more subtle to determine the parameter Ug. In Appendix |B] (around Eq. (IB.16P ). the 
following was shown. Let us we choose the tortoise coordinate to be related to r near the 
horizon as 

- ^ logf^^V (6.12) 



AtiT "\Lh, 
where is defined through the following integral 

for r Tij. Then the solution /aj(r) to the wave equation (I2.42p . satisfying a normalizable 
boundary condition at infinity, will have the form 

L{r) ~ e*^"* - e-'"^^' (6.14) 
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for small uj. More precisely, we have 



Ur) ~ e^"^* - e^-"e-^"^% = 0{uj^). (6.15) 

Now, let us we impose some boundary condition at r^, such as a Neumann boundary condition 
dr^fui = 0, then the frequency u gets discretized in units of Au = vr/|rj|. Note that, if 
a^j = 0{uj) as a; — 7- 0, then the coefficient of the 0{u) term will affect the value of Au; this is 
why f l6.15p was important. This motivates the following choice for the minimum frequency: 

Us ~ Au ~ ~ --J r- ~ -: - . (6.16) 



/31og(^) /31og (^^htirH)hrirH 



Substituting in the above expressions for s, cus, the mean-free-path time (16. 7p is 

^mfp ~ — — 7-7^, V = k^\/ht{rH)hr{rH) , (6.17) 

Tlog(?7VA) ^ 

where A is the '"t Hooft coupling" defined in (15.131) . Note that the nontrivial effect of charge 
only enters through the logarithm and hence the dependence of tmfp on it is very mild in the 
strongly coupled case A ^ 1. 



6.2 Application: STU black hole 

The AdS/CFT correspondence has been successfully used to extract the properties of field 
theory plasmas. A particularly interesting case is a 4-dimensional charged plasma, because 
it is relevant for the experimentally generated quark-gluon plasma with net baryon charge. 
One notable situation to realize 4-dimensional charged plasmas in the AdS/CFT setup is the 
spinning D3-brane, which in the decoupling limit gives c? = 4, A/" = 4 SYM with nonvanishing 
i?-charges. We can have three different i?-charges corresponding three Cartan generators 
of the SU{4) = S0{6) i?-symmetry group. As already mentioned in subsection 12.31 on the 
gravity side this corresponds to a Kerr black hole in AdSs x with three angular momenta in 
the directions [3ZIEH]- Upon compactifying on S^, this reduces to the so-called STU black 
hole of the five-dimensional supergravity [35l|36]. From this five- dimensional perspective, the 
STU black hole is a non-rotating black hole with three U{1) charges. There has been much 
study pTl[39| H0 [l47f[50] on the properties of the i?-charged field theory plasma using the STU 
black hole. Here, we would like to apply the machineries we have developed in the previous 
sections to the computation of the mean-free-path time for the Brownian particle in i?-charged 
plasma dual to the STU black hole. 
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6.2.1 The STU black hole 

The 10-dimensional metric of the STU black hole is given by [35] 



dsl 



fir) 



dfif + fif ( dipi + — 



^dsi + ^yx-' 

-de + v}/' {^ + T\dx'f^ 



(6.18) 



m 



H = H1H2H2 



Xi = H;^H^'\ A' = J- (1 - H-^)dt, A = Y^X.iil 

V 9'i 



/ii = sin 6'i , /i2 = cos 61 sin ^2 , /^s = cos Oi cos 6^2 

with i = 1,2,3. Here, X^ , / = 1,2,3 are spatial directions along the boundary and / is the 
AdS radius. The four parameters m, qi are related to the mass and three electric charges of 
the STU black hole. It is convenient to introduce the dimensionless quantities 



— „ , i — 1,2,3. 



(6.19) 



The horizon is at r = th where th is the largest solution to /(r) = 0. The latter equation 
relates m to rn and 



4 4 

m = "^nirn) = ^(1 + + ^2){l + ^3)- 



The Hawking temperature is given by 

Th 2 + Ki + K2 + K3- K1K2K3 



T 



(6.20) 



(6.21) 



From the five-dimensional point of view, the STU black hole is electrically charged under the 
gauge fields and the associated chemical potentials are 



Here k\ = SttG^ is the five- dimensional Newton constant and 



(6.22) 



VsB 7r3/5 



7cf 

2PP'' 



(6.23) 



^°The horizon of the STU black hole can be either S^, M.^, or H^, but we are focusing on the planar 
case, corresponding to a charged plasma in flat M.'^. 
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where N is the rank of the boundary gauge theory. For expressions for other physical quanti- 
ties, such as energy density, entropy density, and charge density, see e.g. [51]. From thermo- 
dynamical stabihty, the parameters restricted to the range [H2] 

2 — — ^2 — + > 0. (6.24) 

We can shift the gauge potential so that its value on the horizon is zero: 

A'{r)=A\r)-A\rH). (6.25) 

If we accordingly shift the angular variable by 

^,=^, + Ai{rH) (6.26) 

then the metric 06.181) becomes 

dsl, = VAdsl + ^ ^Xri[d/.2 + ^2^d^^ + Aynn (6.27) 

6.2.2 Background configuration 

We first want to find a background configuration of a string in the 10 dimensional geometry 
(16.181) or (16.271) . so that we can start expanding the Nambu-Goto action around it. If we 
restrict ourselves to configurations with trivial 6a dependence, the relevant line element can 
be written as 

3 

ds^ = -a df + I3dr'^ + -f{dX^f + ^ ei(#i + (pidtf. (6.28) 

i=l 

Here a, /3, 7, e^, 0, are functions of r and 9a which can be read off from (I6.27p . For example, 
a = A^/^/7{~^/^. Parametrize the worldsheet by t,r and take the following ansatz: 

X^(t,r) = 0, ^i{t,r) = Uit + ipi{r). (6.29) 

The string is straight in the AdSs part of the spacetime. On the other hand, the angular 
momenta in the directions are expected to drag the string in these directions and (jJi,ipi 
correspond to nontrivial drifting/trailing of the string p!T | l39 | H0]. The Euler-Lagrange equa- 
tion for ip{r) states that = OL^g/ d{dr'ipi) = dL-^c/dipi is constant along the string. The 
quantity tt^^. corresponds to the infiow of angular momenta (or, from the five- dimensional 
point of view, electric charges) from the "fiavor D-brane" at the UV cutoff r = Tc, and how 
to choose them depends on the physical situation one would like to consider ^4Q. Here, let 
us focus on the case where the string endpoint on the "flavor D-brane" is free and there is 
no inflow, i.e., vr^ = 0. This corresponds to a boundary Brownian particle neutral under the 
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i?-symmetry. This is physically appropriate because we want to compute the random force 
correlators unbiased by the effects of the charge of the probe itself. It is not difficult to see 
that setting vr^. = leads to = by examining the Euler-Lagrange equations. 

Let us next turn to the angular velocity Wj. Given Lpi = 0, the induced metric on the 
worldsheet is 

3 

ds'^i^a = -adt^ + Pdr"^ + ^ e,(wi + (t)ifdt^. (6.30) 

i=l 

The determinant of this induced metric is 

3 

det^f oc -a + ^ei(a;j + (6.31) 

1=1 

This must be always non-positive for the configuration to physically make sense. This con- 
dition is most stringent at the horizon r = th where a oc / = 0, 0i = A\{rH)/l = 0. So, we 
need 

^ eiOjj < 0. (6.32) 

i 

Since > 0, this means that 

Ui = 0. (6.33) 

Namely, the background configuration is simply 

X^(t,r) = ^,(t,r) = 0. (6.34) 

Note that the angular motion is trivial only in the ipi coordinates and in the original ipi 
coordinates there is non- vanishing angular drift. 

So far we have been treating 6a as constant. However, this is not correct and an arbitrary 
choice of 6a will not satisfy the full equations of motion. Below, we will consider the following 
three cases: 

(i) 1-charge case: Ki = k 7^ 0, ^2 = ^^3 = 0; 61 = 7r/2, 

(ii) 2-charge case: ki = 0, ^2 = = k 7^ 0; 61 = 0, 

(iii) 3-charge case: Ki = K2 = % = k 7^ 0; 61,62: arbitrary. 

It can be shown [ID] that the above values of 6a are necessary for all the equations of motion 
to be satisfied. These values make sense physically since, if the angular momentum around 
an axis is nonvanishing, the string wants to orbit along the circle of the largest possible radius 
around that axis. This is achieved by the above choices of 6a. 
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6.2.3 Friction coefRcient 



Before proceeding to the computation of tlie mean-free-patli time, let us clieck tliat tlie low- 
frequency friction coefficient for the STU black hole that we can compute using the formula 
fl2.48l) is consistent with the result found in the literature [30] • In the present case of the 
metric (I6.28p . the formula (12.481) gives 



li{uS) = \0{uj). 6.35 

On the other hand, the drag force computed in |30] is|ll 

_7(n«^ (6.36) 

where v is the velocity of the quark and r^g is the solution to a(r^s) — f ^7(r^s) = 0. In the 
non-relativistic limit, f — )■ 0, the admittance read off from fl6.36p should become the same as 
the low- frequency result (16.351) . Using the fact that r^g — ru and p = mv in the f — )■ limit, 
it is easy to see that (I6.36P indeed reproduces the admittance (16.351) . 

6.2.4 Mean-free-path time 

For the three cases (i)-(iii) described above, let us use the formula (16.171) and compute tmfp- 
Consider the n-charge case in = 1,2,3). For the background configuration (16.340 . the 10- 
dimensional metric of the STU black hole (I6.27P induces the following metric: 

ds^ = -fH-'^il - f-'H^A^,)de + H'^[^j + r^(dX'f^ , (6,37) 

r-^ + q r jj + q J r"^ 



At = ^ ( - ) , h = i + ^, (6.38) 



where q = nrj^- Here, in addition to the t, r part, we kept the part of the metric (16.340 also, 
because we would like to consider the transverse fiuctuations along directions. Comparing 
this metric with the general expression (I2.36p . we find 



ht = H-^{1- f-^H'Ai), hr = H^, G = r'H^. (6.39) 
Therefore, from (16.170 . 



^^This is the drag force for the "non-torque string" of [3D] which corresponds to no inflow of at the flavor 
D-brane; see the discussion below (j6.29l) . See Refs. [TTl|3nil3n] for the relation between the strings with and 
without inflow. 
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Figure 4: Behavior of t] versus k, for 1-charge (solid red), 2-charge (dashed green), and 3- 
charge (dotted blue) cases. The range of k, is determined by the thermodynamical stability 
fl6.24p to be K, < 2 for the 3-charge case while there is no restriction for the 1- and 2-charge 



cases. 



The computation of r], particularly Lh in it, is slightly complicated. So, we delegate the 
details of the calculation to Appendix |E] and simply present the results. For 1-, 2-, and 
3-charge cases, rj is given respectively by 



V 



Air 

2 + K 

2n 



exp I -2(2 + k) f 

I J oo 

-4Vl + K [ 

J oo 




- 1 

dp 



exp{ -2(l + /t)2/2(2 



'2 ^ l)(p2 + 1 + 2k) 
^ dp 



(6.41) 
(6.42) 



X 



^(p2 + 1 + ^ _ ^2)(p4 + (1 + 3^)p2 _ ^3) (1 + fi:)3/2 (2 _ «:) 



(6.43) 



The small k expansion of rj is presented in flE.5p - flE.7p . 

In Figure HI we have plotted the behavior of rj as we change k. Because rj appears in the 
denominator of the expression for tmfp, "we observe the following: for the 1- and 2-charge cases, 
^mfp gets longer as we increase the chemical potential keeping T fixed, while for the 3-charge 
case, tmfp gets shorter as we increase the chemical potential keeping T fixed. 

One may find it counter-intuitive that tmfp increases as we increase chemical potential with 
T fixed in the 1- and 2-charge cases, based on the intuition that a larger chemical chemical 
potential means higher charge density and thus more constituents to obstruct the motion of 
the Brownian particle. However, such intuition is not correct. What we know instead is that, 
if we increase the charge with the mass fixed, then the entropy decreases, as one can see from 
the entropy formula for charged black holes. So, if we interpret entropy as the number of 
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"active" degrees of freedom which can obstruct the motion of the Brownian particle, then 
this suggests that tmfp should increase as we increase the charge with the mass fixed. We did 
numerically check that this is indeed true for all the 1-, 2- and 3-charge cases. 

7 Discussion 

We studied Brownian motion in the AdS/CFT setup and computed the time scales char- 
acterizing the interaction between the Brownian particle and the CFT plasma, such as the 
mean-free-path time tmfp; by relating them to the 2- and 4-point functions of random force. 
We found that there is an IR divergence in the computation of tmfp which we regularized 
by introducing an IR cutoff near the horizon. Here let us discuss the issues involved in the 
procedure and the implication of the result. 

The first question that arises is: tmfp is the mean-free-path time for what particle? First 
of all, one can wonder whether this is really a mean-free-path time in the first place, because 
the nontrivial 4-point function was obtained by expanding the Nambu-Goto action to the 
next leading order, which is a relativistic correction to the motion of the bulk string. So, isn't 
this a relativistic correction to the kinetic term in the Langevin equation, not to the random 
force? However, recall the "cloud" picture of the Brownian particle mentioned before; the 
very massive quark we inserted is dressed with a cloud of polarized plasma constituents. The 
position of the quark corresponds to the boundary endpoint of the bulk string, while the 
cloud degrees of freedom correspond to the fluctuation modes of the bulk string. So, we are 
incorporating relativistic corrections to these cloud degrees freedom (fluctuations) but not to 
the quark which gets very heavy in the large m limit and thus remains non-relativistic. 

So, what is happening is the following. First, the constituents of the background plasma 
kick the cloud degrees of freedom randomly and, consequently, those cloud degrees of freedom 
undergo random motion, to which we have incorporated relativistic corrections. Then these 
cloud degrees of freedom, in turn, kick the quark, which is recorded as the random force F felt 
by the quark. F is non-Gaussian, or has a nontrivial 4-point function, because the cloud that 
is interacting with the quark is relativistic. The quark's motion, which is what is observed 
in experiments, is certainly governed by the non-Gaussian random F and the frequency of 
collision events is given by 1/tmfp- However, it is worth emphasizing that this tmfp is not a 
mean-free-path time for the plasma constituents themselveso 

We focused on the fluctuations in the noncompact AdS directions. However, for example, 
in the case of STU black holes, the full spacetime is AdSs x 5*^ and the string can fluctuate in 
the internal directions as well. Let us denote the fluctuations of the string in the internal 
directions by Y, while the fluctuations in the AdS directions continue to be denoted by X. 

^^Ref. 53; estimates the mean-free-path of the plasma constituents to be Zmtp ^ 1/r from the parameters 
of the hydrodynamics that one can read off from the bulk gravity. 



36 



One may wonder if the computations of the random force correlators such as {OxOx^x^x) 
are affected by the Y fields. Here, we denoted the force by Ox to remind ourselves that 
the force is an operator conjugate to the bulk field X. The Y fields do contribute to such 
quantities, because the Nambu-Goto action expanded up to quartic order involves terms of 
the form X'^Y'^ . However, as long as we are interested in quantities with all external lines 
being Ox-, such as xO xO xO x) ■, they only make loop contributions, which are down by 
factors of ol . Therefore, our leading order results do not change. 

In the present paper, we focused on the case where the plasma has no net momentum. 
More generally, one can consider the case where the plasma carries net amount of momentum 
and insert a quark in it. The Brownian motion in such situations were studied in 
(see also [2^) in AdS/CFT setups. It is interesting to generalize our computation of ttmp 
to such cases. Note the following, however: in general, in the presence of a net background 
momentum, the string will "trail back" because it is pushed by the flow. Unless one applies 
an external force, the string will start to move and ultimately attain the same velocity as 
the background plasma. This final state is simply a boost of the static situation studied in 
the present paper. So, the result of the current paper applies to this last situation too (after 
rescaling due to Lorentz contraction). 

The resulting expression for the mean- free-path time, e.g. f l5.12p . is quite interesting be- 
cause of the logarithm. As mentioned around (15.141) . this means that the Brownian particle 
is experiencing ~ log A collision events at the same time. Because X (x N, this is reminiscent 
of the fast scrambler proposal which claims that, in theories that have gravity dual, 

~ log degrees of freedom are in interaction with each other simultaneously. 

In our previous paper pT], we claimed that tmfp ~ 1/T based on dimensional analysis, but 
(15.121) says that there is an extra factor which cannot be deduced on dimensional grounds. Of 
course, we have to note the fact that tmfp we computed in the present paper is not the time 
scale of the constituents but of the Brownian particle (see also footnote [12]) . In our previous 
paper [21], we had t°efax ~ m/{T'^\/X), ~ 1/{T\/X) instead, which were nice because if we 
set m ^ T in f^lf^^ we get f^^. In the (I5.12p . this is no longer the case, but now the relation 
between treiax and t^fp is not so simple as we can see from the fact that there is a nontrivial 
A dependence in the impact per collision, P{ijj = 0) (Eq. (I5.12p ). It would be interesting to 
find an improved microscopic toy model which can relate treiax and tmfp- 

Probably the most controversial issue in our computations is the IR cutoff. When regu- 
lating integrals such as (15. ip . we cut off the geometry at a proper distance s ~ away from 
the horizon, assuming that the contribution from the rest of the integral is of the same order. 
This seems physically reasonable, but we do not have a proof. One could also have tried 
to put a cutoff at the point where the backreaction of the fundamental string on the black 
hole geometry becomes important. Since the interaction of the string with the background 
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is suppressed by additional powers of the string coupling constant, the resulting cutoff is 
presumably closer to the Planck length than the string length. 

Related to the above statement, it is interesting to note that the mean- free-path at weak 
coupling 

9ymTH^/9ym) 

has a form tantalizingly similar to flS.lOp . In particular, the log in fl7.ip is coming from an 
IR divergence cut off by non-perturbative magnetic effects [56], while the log in (15.101) was 
also coming from an IR divergence that we regularized by introducing an IR cutoff. It would 
be interesting to study whether there is a relation between the weakly and strongly coupled 
descriptions of the IR divergences and the physical interpretation of the IR cutoffs. 
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A Normalizing solutions to the wave equation 

As explained in subsection 12.21 or more generally in subsection 12. 3[ the normalized modes 
{u^} are proportional to of the form fl2.22p : namely, u^(t,r) oc e~*'^*/(^(r). Here, we fix the 
normalization and derive the expansion fl2.23p or more generally f l2.47p . 

The analogue of the Klein-Gordon inner product for functions f{t, r),g{t, r) satisfying the 
equation of motion (12.420 is [21] 

(/, 9h = ^ n^G {fd,g* - d,f g*), (A.l) 

where S is a Cauchy surface in the t, r part of the metric (12.360 . g is the induced metric on 
S and is the future-pointing unit normal to S. 
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We want to normalize using this norm ( lA.ip . In the present case, there is the following 



simplification to this procedure. Near the horizon r ^ rn, the action (12.411) reduces to 

^ j dtdn[{dtXf-{drXy]. (A.2) 



Ana' 



Therefore, in this region and in the tortoise coordinate system, X is just like a massless 
Klein-Gordon scalar in flat space. Correspondingly, the contribution to the norm (lA.ip from 
the near horizon region is 

f dnifd,g*~djg*), (A.3) 

where as S we took the constant t surface. This is the usual Klein-Gordon inner product for 
the theory (1A.2I) . up to overall normalization. Of course, there is a contribution to the inner 
product from regions away from the horizon. However, because the near-horizon region is 
semi-infinite in the tortoise coordinate (recall that r = th corresponds to = — oo), the 
normalization of solutions is completely determined by this region where the inner product is 
simply (lA.3p . This means that the canonically normalized mode expansion is given by 



^^''^^ = \iWH)[^^i^ [/.We-^-*a. + /.(r)*e-*at] . (A.4) 



where /i^(r) behaves near the horizon as 

U{r) ^ e*""* + e''-e-''^''% r^m (n ^ -oo) (A.5) 

with some 9^^ G M. If we can find such /tj(r), then a,a^ satisfy the canonically normalized 
commutation relation (I2.15p . 

B Low energy solutions to the wave equation 

Here, we study the solution to the wave equation (I2.13p . or more generally (I2.42p . satisfying 
an appropriate boundary condition (the Neumann boundary condition f l2.16p or normalizable 
boundary condition at infinity), for very small frequencies u. We see that the solutions become 
trivial plane waves in the limit. 

The general wave equation fl2.42p can be written in the frequency space as 



V rir ^ \ V f^r 



XJr) = 0. (B.ll 



Very close to the horizon, this becomes 

[uj^ + lQn^T\r - rH)dr{{r - rH)dr)] X^(r) = 0. (B.2) 
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This means that the hnearly independent solutions are 



g±^ = exp 



±«^^ log 



AtcT 



(B.3) 



where Li is a length scale which is arbitrary at this point. The ± signs here correspond to 
outgoing and ingoing waves. We are considering the small u limit but, no matter how small 
u is, we can always consider a region very close to the horizon so that log C'^l^ ) ~ ^(■'-)' 
namely < ^-at^t/u) ^ such a region, we cannot expand the exponential and should keep 



the full exponential expression fIB.Sp . In other words, the precise limit we are taking is 

'r — Th 



a; 0, 



: fixed. 



(B.4) 



Now, consider the region not so close to the horizon. For small w, we can ignore the o;^ 
term in (IB.ip . obtaining 



= Bi + B2 



dr' 



f{r')G{r') V ht{r 
where -Bi,i?2 are constant. For r ^ r^, this gives 

B2 , fr-ru 



(B.5) 



Xuj = Bl + 



log 



Here, we defined the constant Lh by 

dr / hr 



L 



H 



+ 0{r - Th) 



(r ~ ru). 



(B.6) 



fG V ht A-KTG^ru) 



L 



H 



log( ]+0{r-rH). 



(B.7) 



Because it will turn out to be convenient to choose Lh = Li, we will set L^ = Li henceforth. 
On the other hand, for large r, (IB.Sp gives (assuming the large r behavior f l2.39p . ( I2.40p of 
functions ht,hr,G), 



(B.8) 



We can determine Bi, B2 by comparing these small-frequency solutions between the very- 
near-horizon region and the not-so-near-horizon region. For small frequencies u, (1B.3P be- 
comes 



47rT 



log 



(B.9) 



Comparing this with (]B.6|) . we determine 



B^ = l + 0{u^), B2 = ±tuG{rH) + 0{u^). 



(B.IO) 
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Therefore, the hnearly independent (outgoing/ingoing) solutions are 



exp 



1 ± 



±«- — — log 



A-kT 

iujG{rH] 
3r3 



(B.ll) 



The general solution is given by the linear combination of the outgoing and ingoing 
solutions g±i^. If we want to construct a normalizable solution that vanishes as r — )■ oo then, 



from the r ^ th behavior of fIB.lip . the linear combination to take is 

li^(norm) _ 



(B.12) 



If we did not take Li = Lh, the two terms would be multiplied by exp[=Fi^^ \og{^)\ respec- 
tively. Note that our expressions are correct up to 0{uj'^) terms. The near-horizon behavior 
of this is 



00 ( r — Th 

l- — =: log 



AttT 



L 



H 



exp 



00 r — Th 

-I- — =: log 



AtiT 



(B.13) 



Therefore, if define the tortoise coordinate to have the following behavior the horizon: 



1 



lot 



r — Th 

then the near-horizon behavior flB.lSp simply becomes 



j^(norm) ^ ^jur, _ ^- 



(B.14) 



(B.15) 



Let us elaborate on this point slightly more. In the near horizon region, in general we can 
have 



ji^(norm) ^ ^iu 



(B.16) 



where a^^ is some phase. The fact that (IB.lSp is correct up to O{oo) means is that, if we take 
r* to be given by fIB.Mp . then = O{oo^) as a; — )■ 0. In particular, unless we choose Lh to 
be the one given by (]B.7p . the a; — ?■ behavior of a^j will contain an O{oo) term. 

Next, let us consider imposing a Neumann boundary condition drX = at r = ^ r/^ 
instead. Set the general solution to be 



then the Neumann boundary condition X^(rc) = gives 



(B.17) 



C 



O{oo^ 



1 + 0{uj) 



(B.18) 
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where we used the second equation in (IB. 111) . Comparing this result with (12.211) and fl2.24p . 
we find that the modes f^i satisfying the Neumann boundary condition are given by, at low 
frequencies, 



9M) = e*"'% 0^ = 0, Ur) = e*"'^* + e"^'^''*. (B.19) 

le 

S 



This is consistent with the explicit result for AdSa in fl2.30p . fl2.32p . So, for very small u, t 
solution fuj{r) is a simple sum of outgoing and ingoing waves, which are just plane waves 
Because g±{r) — )■ 1 as r — )• oo, we have 



Ur = r,)^2. (B.20) 



Because of the 0{u) ambiguity in fIB.lSp . 9^^ = 0{u) as a; — )■ (cf. comments below flB.16P ) 



In the tortoise coordinate, the limit ( ]B.4p we are taking can be written as 

u ^0, /3, un ■■ fixed. (B.21) 

C Various propagators and their low frequency limit 

The quadratic action for a string embedded in the AdS^ black hole spacetime 

ds^ = -hAr)fir)de + ^^dr^ + G{r)dX\ (C.l) 



(Eq. ([23SD) is given by 



S'o = / dtdr 

Ana' 



(C.2) 



We would like to regard this system as a thermal field theory at temperature T, and derive the 
relation among various propagators (Green functions) and the solutions to the wave equation. 
We present the result for the general metric (12. 36 p . but if one wants the results for the simpler 
neutral case (I2.10p . set / = r'^h, ht = hr = 1, G = r^. 

Let us define Wightman, Feynman, retarded, and advanced propagators as 

Dw{t-t',r,r') = {X{t,r)X{t',r')), 

D^it-t',r,r') = {r[Xit,r)Xit',r')]), 
D,,,it-t',r,r') = eit-t'){[Xit,r),Xit',r')]), 
DMt - t', r, r') = -e{t' - t) ([X(t, r), X(t', r')]). 



13 



For related observations on the triviality of the solution in the low frequency limit, see [5 7) . 
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We impose a Neumann boundary condition for X{r,t) at r = r^, so the propagators satisfy 
the same Neumann boundary condition. Using the wave equation 



and the canonical commutation relation 




Gfdr 



X = 



(C.4) 



[X{t, t) , dtX{t, r')] = 2ma'J^L 5{t - r') 



(C.5) 



we can show that these propagators satisfy 



G 



htf y/Jhh. 



zdr 




Gfdr 



Dwit-t',r,r') = 0, 



G o2 



^F,ret,adv(i - i', r') 



27cia' 



(C.6) 

5{t-t')5{r-r'). (C.7) 



where y^—g = \fhtK.. 

As in flA.4l) . the field X can be expanded as 



X(t,r) 



2'Ka' 



duj 1 
2tx J2Lj 



[/.(r)e-''"*a^ + /.(r)*e^'^*ay , 



(C. 



where 



and guj^r) behaves near the horizon as 



(r ~ rH) 



(C.9) 



(C.IO) 



The phase Qui is determined by the Neumann boundary condition aX r = Vc that f^j satisfies. 
Since the system is at temperature T, the expectation value of a, a"'' is given by fl2.25p . It is 
then easy to show that the Wightman propagator can be written as 



Dwiuj,r, r') 



G{rH) 2w(l - e- 



(C.ll) 



where = f*. 

We would like to express other propagators -Dadv,ret,F in terms of fuj.gui- Note that 



Df{uj, r, r') = Dw{(^, r, r') + L'adv(w, r, r') = r', r) + L'ret(w, r, r'). (C.12) 
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Because we have already obtained Dw in (IC.lip . if we know one of Dp, Dj-^t, and -Dadv, we 
can obtain all other propagators. Here, let us consider -Dadv From flC.7p . Dadv(i^, r, r') should 
satisfy 



-OJ H . Or 




Gfdr 



-Dadv(w,r, r') 



2'Kia' 



-9 



6{r — r'). 



(C.13) 



htf yhtK \ V K 
If r ^ r', this is the same as the wave equation that satisfy. Therefore, take the ansatz 

Dadv(u;, r,r') = A[e{r - r')g^{r')Ur) + e{r' - r)g^{r)Ur% (C.14) 

This satisfies the correct boundary condition (Neumann) at r, r' = and furthermore satisfies 
the purely outgoing boundary condition at the horizon, which is appropriate for an advanced 
correlator. Using the fact that both /, g satisfy the wave equation, we find 



hi VhtK 
Therefore, 

A = 27ria- 




hr 



, jlVr 1 



Gfdr 



Dadv = ^^^'^ '''^ .f^Gf {g^drU - Org. U) ■ (C.15) 



2TTia' 



(C.16) 



ht Gf g^jdrfu - drguj fu Girn) {gujdrju - dr^g^ fio) ' 

Using the wave equation for fu},gui, it is easy to show that this expression does not depend 
on r. By taking r ^ and using (lC.9p . (IC.lOp . 



A 



So, the advanced propagator is given by 



Tia e 



G{rH)co 



(C.17) 



na'e 



[e{r - r')g^{r')f^{r) + e{r' - r)g^{r)f^{r')] . 



(C.18) 



CM a; 

In the low frequency limit, the expressions for the propagators simplify, as we saw in 
Appendix [HI The precise limit we are considering is f lB.2ip . First, the Wightman propagator 
(IC.lip becomes, because of (IB.lQp . 



Similarly, the advanced propagator (IC.lSp becomes 



(small w) 



(C.19) 



Dadv(uj,r, r') 



GirH)uj 



Girn) 



(small to). (C.20) 
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Using the relation (1C.12|) . the Feynman propagator is 



DF{uj,r, r') 



G{rH) 



(small u). 
(C.21) 

In particular, consider the case where one of the points is at the UV cutoff, r' = r^. From 
flB:20|) . we have 



Dw{uJ,r, r-c) 



27ia' 



(C.22) 



D Holographic renormalization and Lorentzian AdS / CFT 

In this Appendix, we discuss how to compute correlation function using the AdS/CFT dic- 
tionary for the total force F which is dual to the worldsheet field X. As we explained in 



subsection \A.2\ this involves holographic renormalization (see e.g. [2^) of the worldsheet ac- 
tion. Furthermore, if we want to compute real time correlation functions in a black hole (finite 
temperature) geometry, we should apply the rules of Lorentzian AdS/CFT [271128]. 

D.l Holographic renormalization 

First, let us consider the holographic renormalization of the worldsheet action. For this, 
only the asymptotic behavior of the action near the boundary is relevant. Therefore, as the 
background geometry, we can consider the Poincare AdS geometry obtained by setting T = 



ds' 



-rdr + ^ + r\dX 



I\2 



for which the worldsheet action becomes 

'S'barc = 'S'o + Sint, 
1 



dtdr{X' -r^X'^) 



Si 



1 



int 



dtdr{X' -r^X 



4 vi2\2 



(D.l) 

(D.2) 
(D.3) 



Ana' Jy. IQna' j^, 

Here, we considered only one of the polarizations, say X^, and denoted it by X . E is the 
worldsheet, 



S = {{t,r)\te M, < r < rj. 



(D.4) 



and ' = dt, ' = dr- is a UV cutoff. For computational convenience, let us we rescale 
X — )■ y/2'n'a' X and set n = ira', so that 
1 



So 



4 v/2\ 



dtdriX' -r'^X 



Sir 



4 



dtdr{X^-r^X 



4 v/2\2 



(D.5) 
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The equation of motion is 

-d^X + dr{r%X) = K[-dtiHdtX) + driHr%X)], H = -X^ + r^X'\ (D.6) 
Let us solve the equation of motion (1D.6P by expanding X{t,r) in the couphng k as 



X{t, r) = Y{t, r) + KZ{t, r) + 0{k^) 
and furthermore expanding Y, Z around r = oo as 

Henceforth, we will ignore quantities of The expansion for X itself is 



(D.7) 



(D.8) 



A [t, r) = X(o){t) H \ \ h 



= Uii) + KZ{i). 



(D.9) 



By substituting this expansion into flD.6p and comparing coefficients, one readily finds that 
the following is a solution: 



2/(0) = any = J, 
Z(o) = 0, 
The expression for X is 

J, = 



2/(1) = 0, 

0, 2^(2) 



vm = 



-¥ 



2/(3) = any, 
Z{3) = any. 



X 



(0) 



Xi2) 



--J-kJ^J, 
2 



X 



(3) = any. 



(D.lOa) 
(D.lOb) 

(D.ll) 



Note that X{r,t) — t- J{t) as r — )■ oo; namely, J is the non-normalizable mode which can be 
thought of as a source for the dual operator Ox = F on the boundary. On the other hand, X(3) 
is the normalizable mode which roughly corresponds to the expectation value of the operator 
F. We will make this latter statement more precise below. 

If we plug the solution f lD.lip into the action ( ID.Sp . we obtain the following on-shell action: 

1 



5 



bare, on-shell 



d^xJJ.J + 



dt 



-rJJ — nrJJ J 



K 
—1 

4 



-rjrj 



dt 



as 



IrJJ-^rjj'j 
2 4 



dt 



as 



+ (finite) 
+ (finite). 



(D.12) 



In going to the second line we performed the r integration, and in going to the last line we 
integrated by parts. This is divergent, but the divergence can be canceled by introducing the 
following counter terms: 



S. 



ct 



dt a/— 7 



as 



(D.13) 
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where = —^d'f is the Laplacian for the metric 7 induced on the boundary r = r^. Likewise, 
(V-yXy = — ^((9(X)^. If we define the metric 7' induced on the boundary of the worldsheet 
at r, then 7^^ = — r^(l — X^) and J dt a/— 7^^ reproduces flD.lSp (also recall that we have 
rescaled X — )■ y/2iTa' X). 

To remove the divergence from the "bare" action fID.Sp . we take Sren = Share + Set as our 



total action. The on-shell variation of this total action evaluates to 

55,en,on-shell = j ^^dt,/^ (^-T^d^X + D^X) + KT^ [{yXfd^X + ?>{y ^XfU^^X]^ 5X , 

(D.14) 

where dn is the normal derivative with respect to the worldsheet boundary dH. Therefore, 

5'ren , on-shell 



5J 



^(^G{dnX + U^X) + kG"^ [{VXfdnX + 3(V^X)2n^X] ^ 
3x(3)(l + /tj') + C(l/r). (D.15) 



In the second equality, we plugged in the explicit expansion (ID. lip . Therefore, by the GKPW 
rule OIH], the expectation value of the operator Ox = F dual to X in the presence of source 
X(o) = J is given by, up to 0{k'^) terms. 



{F)j = 3x(3)(l + kJ') = 3|/(3) + 3/€ (^2;(3) + 1/(3) J'j . (D.16) 

The term may appear strange, but we will see that this term gets canceled in the final 
expression for the 4-point function. Actually, there is a further contribution to f lD.16p . but 
we will discuss it later (see below f lD.40p ). 

Although our discussion above was based on the pure AdS space (ID.ip for the simplicity 
of the argument, the final expression flD.16p is valid for general asymptotically AdS space, in- 



cluding the AdS black hole fl2.10p . Below, we will use f lD.16p to compute correlation functions 
for the AdS black hole background (12.101) . 

D.2 Propagators and correlators 



To compute the expectation value {F) j using the formula f lD.16p . we need to know X(3) = 
l/(3)+fi;Z(3)+(9(K^). This can be determined if we know the propagators that satisfy appropriate 
boundary conditions in the inside of the AdS space as we discuss below. 

If we substitute the expansion (1D.7P into the wave equation (12.131) and compare the coef- 
ficients, we obtain 

[-h'^d^t + dr{r%dr)]Y = 0, (D.lTa) 
[-h-^d^ + dr{r^hdr)]Z = p, (D.17b) 
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where we are now considering the AdS black hole spacetime fl2.10p and the "source" p is 
defined by 



P 



-OtiHoh-^dtY) + driHor^hdr-Y), Hq = -h-\dtYf + r^/i(5,F)^ 



(D.18) 



We solve (]D.17a|) under the asymptotic condition Y{r,t) J{t) as r — )• oo and (]D.17b|) 
under the condition Z{r,t) — )■ as r — )• oo. Let us solve these using propagators. First, 
let K{r,t\t') be the boundary-bulk propagator, namely the solution to the zeroth-order wave 
equation (]D.17a|) satisfying the boundary condition 



K{r,t\t') ^ 5{t -t') as r ^ oo. 
Then the solution to (ID.17aP is 



F(t,r) 



dt'K{r,t\t')J{t'). 



(D.19) 



(D.20) 



From this, we can read ofi' ?/(3) as 



|/(3)(t)= / dt'[K{r,t\t')]\,.-zJ{t') 



(D.2r 



where [ j^-s means to take the coefficient of the term in the 1/r expansion. 

Let us move on to the next order equation (lD.17b|) to determine Z(3). Let D(r, t|r',t') be 
the bulk propagator, namely the solution to 



[-h-^dt + dr{r^hdr)\D{t, r\t', r') = 6{t - t')6{r - r') 



(D.22) 



that vanishes as r, r' — t- oo. Then the solution to the next order equation flD.17bl) can be 
written as 



Z{t,r) 



dt' dr' D{t, r\t', r')p{t , r ) 



(D.23) 



It is easy to see that the Z given by flD.23p has the expected behavior fID.lObI) . To see it, 
let us explicitly construct the bulk propagator satisfying (1D.22|) . or in the frequency space, 

[h-^uo"^ + dr{r^hdr)]D{uj, r, r') = 6{r - r'). (D.24) 

The solution to this can be constructecj^ from the solution to flD.17al) . which can be written 
in the frequency space as 



[h-^io' + dr{r^hdr)]Y^ = 0. 



(D.25) 



^■^The following argument is analogous to the one given around (jC.13p . 
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As discussed above the equation (12.181) . this wave equation (1D.25P has two solutions; let us 



denote them by 0±aj('")o These are related to each other by 0;^(r)* = (t>-uj{r). As one can 
see from fID.lOap . we can take them to have the following large r expansion: 



0±.(r) = l + ^ + ^ + ..-, (D.26) 
where c±oj are some constants (c* = C-^j). For example, in the AdSs case {d = 3), 



For r 7^ r', the equation (ID.24P is the same as flD.251) and therefore D{u,r,r') is given by a 



linear combination of 0(^(r) and (f)-^{r). Taking into account the r -r^ r' symmetry, the bulk 
propagator D can be written as 



€ir)<l>^iryir - r') + C(r')0<(r)^^(r' - r) . (D.28 



(D.29) 



Here A is constant and we defined 

0>(r) = Ur) - <P-Ur) = ^^^^^^ + 0{r-'), 

4><{r) ^ aUr) + (1 - «)0-.(r) = 1 + |^ + ^ ^^^3" ""^'"^ + ^(^"')- 

The fact that 0^(r) — t- as r — )■ correctly gives the asymptotic condition for D, namely 
D — i- as r, r' — )■ 00. On the other hand, we do not specify the boundary condition of D as 
r, r' — )■ vh- The unknown number a parametrizes possible boundary conditions which is to 
be determined by some physical requirement. But we leave a arbitrary and therefore (]D.28|) 
is valid regardless of the boundary condition. Because — )■ 1 as r — > 00, it is actually 
equal to the bulk-boundary propagator in the frequency space; 

<P<{r) = K{co,r). (D.30) 

By substituting flD.28p into the equation flD.24p . we obtain 

^=r4M0>(9.0<)-(9.0>)0<] ^^-^^^ 

(this is the same as (IC.16P ). Since this does not depend on r (see below (1C.16P ). by taking 
r — 00 and using the asymptotic behavior flD.29p . we find A = (c^ — c_^)~^. Therefore, the 
bulk propagator is found to be 

Diu,r,r') = (c^-c_^)-i[0>(r)0<(r')e(r-rO + 0>(r><(r)^^(r'-r)], (D.32) 



^^(f>±uj{r) are different from g±ij{r) defined around (|2.18p only by normalization; (('±^{1") 1 as r — )■ 00, 
while g±uj{r) — > e^^'^^' as r — >■ rc (r* —00). These agree in the small ui limit. 
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where we used (lD.30p . The r — )■ oo behavior of this is, using the asymptotic behavior (lD.29p . 
D(u,r,r') = -^K(u,r')e(r-r')-^e(r' -r) + 0(r-'') (r ^ oo). (D.33) 



r')9{r - r') - -^9ir' - r) + 0(r-^) (r ^ 



Using fID.lOap . we can show that the source p (defined in Eq. flD.lSp ) goes as p = 2.PJ + 
0{r^'^). Then, from (]D.23p and ( 1D.33P we can read off ^(3) as follows: 



Z(3)(t) = lim 



dt'dr'K{t',r'\t)p{t',r') + -rj{tfj{t) 



(D.34) 



The second term cancels the divergent contribution corresponding to Z(2) in flD-lObp . 

So, we succeeded in expressing 7/(3), ^(3) appearing in the formula (1D.16P using propagators; 
the resulting expressions are flD.2ip and (ID.34p . Using these, we can compute the boundary 
correlators for F. First, at the first order in n that we are working in, the 2-point function 
gets contribution only from 7/(3) in (ID.2ip and 



{T[F{t{)F{t,)]) 



In the frequency space. 



5 



6J{t2 



S 



.7=0 



5J{t2 



-y(3)(ti) = 3[ir(ti,r|t2)]|.-3. 



{F{uji)F{lo2)) = 27r5(wi + L02) 3K(w2,r)|^-3. 



(D.35) 



(D.36) 



To obtain 4-point functions, we take functional derivatives of (1D.16P three times. There- 
fore, only the second term 3k(2;(3) + y(^yP) in (ID.16P is relevant for the computation. Let us 
write the source p appearing in f lD.34p as 



p = dtp"" + drp\ p* = -Hoh-^dtY, p'' = H^r^hdrY. 
Then, by partial integration, f lD.34p becomes 



(D.37) 



Z(3)(t) = lim 



dt'dr' 



dt' 



p\t\ r')dt>K{t', r'\t) + p'{t', r')dr>K{t', r'\t) 



K{t',r\t)p^{t',r) - K{t',rH\t)p'{t',rH)\ + -rj{tyj{t)}. (D.38) 



We dropped the boundary terms at t = ±00. The first term in the second line can be 
evaluated using the expansion 

K{t',r\t) = 6it-t') + Oir-^), p'{t,r) = rj^J + Sy^^^j^ + Oir-^). (D.39) 

As a result, in the combination appearing in flD.16p . the term involving y{z)J'^ cancels out: 



K lim 



dt'dr' 



p\t', r')d,K{t', r'\t) + p\t', r')dr>K{t', r'\t) 



dt'K{t',rH\t)p'{t',rH) + rdt[j{tY 



(D.40) 
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The second last term in (1D.40P gets canceled by the extra contribution alluded to below 
flD.161) . Let us now discuss what this extra contribution is. The on-shell variation of the 
action, which we used to compute the expectation value {F)j, is given by flD.14p . Because 
we are regarding the region r^^ < r < as our spacetime, there actually is contribution 
from the "boundary" r = rn to this expression. In the AdS black hole spacetime, this extra 
contribution to 6Sren on-sheii becomes 



SSr. 



shell ^ - / dtr%{drY + KHQdrY)5Y, 

Jr=rH 



(D.41) 



where we dropped 0{k'^) terms and "d" means that the left hand side includes the expression 
on the right hand side. Note that, because the counter term Set (EH]) was added only for 



the boundary at infinity, the second and the fourth terms in flD.14|) did not contribute to this 
expression. Since /i — )■ as r — r^, this becomes 



SSr, 



-shell D — K 



dtr^hHodrY 5Y 



(D.42) 



r=rH 



(note that Hq involves h ^). Therefore, by taking functional derivative, we find that there is 
the following extra contribution to {F)j: 

'^'S'l-en, on-shell 



6J{t) 



D-Kl dt'r'^hHodrY K{t',r'\t). 

' r'=ru 



(D.43) 



This precisely cancels the second last term in (ID.40p . Therefore, the terms relevant for 
computing 4-point functions is 



rH 



{F{t))j DK\im\ / dt'dr' jt{t' , r')dt'K{t', r'\t) + >(t', r')dr'K{t', r'\t) + rdt[J{t) 



(D.44) 



By taking functional derivatives of flD.44|) . we find that 



5^ 



K 



4^ 



dt dr 



6J{t2)6J{h)6J{U 



J=0 



h 



(D.45) 



where the r — oo limit is understood. Also, Ki = K{t,r\ti) and the summation is over 
permutations (ijkl) of (1234). The expression in the Fourier space is 



[ i^i^jK.K, + r%K[K';) {cOkUiKkKi + rW^K[) - Qru^uj2UJ3uX (D. 

V. perm J T fj ) 

(ijkl) 



46) 
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Figure 5: The bulk geometry M = Mi + M2 + M3 that "fills in" the boundary contour 
C = Ci + C2 + C3. For d > 3, the Penrose diagrams for the Lorentzian patches drawn above 
are not accurate because the zigzag singularity lines must actually be not horizontal but bent 
inwards l58l. 



where now Ki = K{ui, r). Note that the first term in f lD.46p is the expression for the 4-point 
function we would obtain from the naive GKPW rule. The last term is there to cancel the 
UV divergence coming from the first term due to the fact that Ki = 1 + (9(r~^). 

D.3 Lorentzian AdS/CFT 

So far we have not fully taken into account the fact that our spacetime is a Lorentzian 
spacetime, for which we have to use the Lorentzian AdS/CFT prescription [27|l28]. 

On the boundary side, to compute real time correlators, we have to take the time to run 
along the contour on the complex place, as we discussed in subsection 14. Ij see Figure |3]on page 
[20I The Lorentzian AdS / CFT prescription is simply to consider a bulk spacetime which "fills 
in" this contour. Then the bulk spacetime will have no boundary and there is no ambiguity 
in boundary conditions (although we have to impose certain gluing condition for fields across 
different patches). Following [28j, we take the bulk spacetime to be the union of three patches 
Mi with i = 1,2, 3, each of which fills in the corresponding contour Ci in (14. 5p . First, we take 
Ml to be the —L <t<L,rH<r<oo part of the Lorentzian AdS black hole (I2.10p . M2 is 
taken to be the same as Mi metric- wise, but the orientation is taken to be opposite to Mi, 
corresponding to the fact that Ci and C2 has opposite orientations. M3 is taken to be the 
Euclidean version of the black hole (I2.10p . 

2 72 
dsl = ^ [hir)dT'' + {dX'Y] + -^dr^ (D.47) 

The Euclidean time r is taken to be < r < /3 where /3 is the inverse Hawking temperature in 
(12.111) . For a schematic explanation of the patches Mi_2,3, see Figure [51 The way that three 
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Figure 6: How to patch together the bulk patches Mi, M2, M3. 

patches Mi^2,3 are glued together is simply the bulk extension of the way that the contours 
^1,2,3 are glued together; see Figure El 

Because now our spacetime is not just Mi but M = Mi + M2 + M3, the action have 
contributions from all of Mi^2,3, just as the boundary f l4.8p . Therefore, the bulk integration 
appearing e.g. in ( ]D.46|) should be now over all Mi, with the signs correctly taken into account: 



G^{uJi, UJ2, ^3, W4) = 2'nK5{uji + + ^3 + W4) 

X 



4 

perm r/f 



— QruiU2U}zU}i J. (D.48) 

Here, i^fabji = K]^ab]{^iif) and i^[ab](w,r) is the boundary-bulk propagator from the boundary 
dMh to the bulk Ma- The second line corresponds to the integration over Mi and the third 
line to the integration over M2. Because we are taking the L — )■ 00 limit, the contribution 
from M3 has been dropped. The counter term —QrujiUj2Uj-iUJi is added only for Mi, because 
the source is inserted only on dMi (_ft'[2i](co', r) vanishes as r — )■ 00). 

Because the spacetime M = Mi + M2 + M3 has no boundary inside, the boundary-bulk 
propagator can be determined without having to worry about boundary conditions. Carefully 
matching the values across different patches following |2ZII2B], we find the boundary-bulk 
propagators as follows: 

= ^^^[-0.(r) + 0-.(r)], (D.49) 
ir[3i](w,r) = — -[-0^(r) + 0_^(r)]. 



where (l)±u){f) is the solution to the wave equation (ID.25P satisfying the boundary condition 



(ID.26p . By substituting these propagators into flD.48p . we can finally obtain the 4-point 
function for F . 



53 



D.4 Low frequency correlators 

We are interested in the low frequency behavior of the correlation functions. As we discussed 
in Appendix [HI the solution 0±t^(r) simplifies in the low frequency limit a412 

0±.(r)~e±*-^*. (D.50) 

If we apply this to flD.48|) and (lD.49p . we obtain the following low frequency behavior: 



G^{ui,U2,U3,U4) —6{uJi+U2 + U3 + U4) 2j i^i + ^j) 

" l<i<7<4 -^-"^ 



h 



i<3t 

+ (higher powers in w), (D.51) 



where we dropped numerical factors. Because we have rescaled X in f lD.Sp . to obtain the 



correlator for F = Ox dual to the original X before rescaling, we have to rescale F — )■ 
Therefore, in the end, the 4-point function for F is 

rO ^2 



1 r r 

G^{ui, U2, 0J3, W4) ~ ^'^(^1 + W2 + + W4) + ^J') / c^r,— e-2^('^»+'^^>*. 

(D.52) 

This is exactly the same as the result fl4.15p that we obtained by a more naive method. 
Namely, this has exactly the same IR divergence as (I4.15p that we studied in section O 

D.5 Retarded 4-point function 

In the above, we computed the time-ordered 4-point function for the force F which turned 
out to be IR divergence. We can also compute the retarded 4-point function using the above 
formalism. As was shown in |2E], for computing retarded correlators, one uses purely ingoing 
boundary condition for the boundary-bulk propagator: 

K,,,{u,r) = (l)_^{r). (D.53) 

If we define 

GUti,h,h,U) = > h >tk> ti){[[[F{U),F{tj)],F{tk)],F{ti)]) (D.54) 



perm 

[ijki) 



then the prescription of [28] gives 

Gf^^{ui, U2, UJ3, UJ4) = 2nK5{uji + U2 + UJ3 + (^i 

X 



1 V r' 

perm ^Tff 
(ijkl) 

— QruJiU2UJ3UJ4 \, (D.55) 



^^Note that the precise hmit we are taking is (jB.2ip . 
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where K^et,i = Kj-^ti^i, r) . The integration effectively becomes only over Mi. 
For definiteness, consider the AdSs case where the retarded correlator is 



r / \r + Th 



For this case, Eq. (]D.55|) gives 



GZ.JuJi,U2, UJ3, U4) = k27i5(uJi H h W4) UJlUJ2UJ3^4: ('^th ^^^^^ — ■ (D.57) 

V Th J 

Note that this is exact; we have not done low frequency approximation. This is both IR and 
UV finite. 



E Computation of 77 for the STU black hole 



In this Appendix, we will compute the mean-free-path time for the STU black hole studied 
in 16.21 The final results have been presented in fl6.40p and fl6.4ip - fl6.43p . 

We will discuss the 1-charge case (ki = k, K2 = K3 = 0) only, because the 2- and 3-charge 
cases are similar. First, the relations fl6.20p . fl6.2ip . and fl6.22p read, in this case. 



m 



T 



(E.i: 



Lh in f l6.40p can be computed as follows. From the definition fl6.13p and fl6.39p for n = 1, 
we obtain 



00 r^f 



1 



1 - n^Al 



dr 



H 



log — h 0[r - Th) 



The integral in the first line diverges as r — )■ rn- We can separate this divergent piece by 
subtracting and adding the term obtained by setting r to th in the square root. Further 
setting p = r/rfj and k, = (? /th-, we have 

dp 



1 y/TT^ 
2r% (2 + k) 



log + 2(2 + /t) 



1 VT 



log 



r + rH 
r — Th 



1 



^{p + 1 + k){{1 + k)p^ + 1) (2 + k) 



+ 0{r-rH) 



(E.3) 



2r|,(2 + /c)"° Lh 

In the second term in the first line, we have set the upper limit of the integral to p — 1 (which 
is equivalent to r — )■ rn) because the integral is now convergent. By comparing both sides, 
we obtain 

1 1 



Lh = 2rH exp | -2(2 + k) j 



v/(p2 + l + K)((l + K)p2 + l) 2 + «: 



(E.4) 
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By using (lE.ip . we obtain the final expression fl6.4ip . For small k, it is easy to expand the 
integrand in f l6.4ip in k, and each integral converges. This leads to the following expansion 
of T] in k: 

(12 - ttW 



r] = e 



-n/2 



2tT — TTK + 



16 



-K 



(E.5) 



This shows that, as n increases with fixed T, the mean-free-path time t^/p increases. 

The 2- and 3-charge cases are similar and we obtain fl6.42p and fl6.43p . The small k 
expansion of i] is 



rj = e 



-tt/2 



-7r/2 



1] = e 

for the 2- and 3-charge cases, respectively 



2n - ^(4 - n)nK + -^vr (52 - 19n + n^) + O U^) 
2 16 

27r + (tt - 3)7rK + ^vr (l40 - 577r + An^) + O {k^) 



(E.6) 
(E.7) 
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